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Abstract 

Topical problems in the physics of and basic facts about 
neutron stars are reviewed. The observational manifesta- 
tions of neutron stars, their core and envelope structure, 
magnetic fields, thermal evolution, and masses and radii 
are briefly discussed, along with the underlying micro- 
physics. 
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1 Introduction 

Neutron stars, the most compact stars in the Universe, 
were given this name because their interior is largely 
composed of neutrons. A neutron star of the typical 
mass M ~ 1-2M0, where M q = 2 x 10^^ g is the solar 
mass, has the radius i? ~ 10—14 km. The mass density p 
in such star is ~ 10^^ g cm~'^, or roughly 3 times normal 
nuclear density (the typical density of a heavy atomic 



nucleus) po = 2.8 x 10 g cm~^. The density p in the 
center of a neutron star can be an order of magnitude 
higher than po- Such matter cannot be obtained under 
laboratory conditions, and its properties and even com- 
position remain to be clarified. There are a variety of 
theoretical models to describe neutron star matter, and 
a choice in favor of one of them in the near future will 
be possible only after an analysis and interpretation of 
relevant observational data using these models. Neutron 
stars exhibit a variety of unique properties (that are dis- 
cussed below) and produce many visible manifestations 
that can be used to verify theoretical models of extreme 
states of matter [1]. Conversely, the progress in theoret- 
ical physics studying matter under extreme conditions 
creates prerequisites for the construction of correct mod- 
els of neutron stars and adequate interpretation of their 
observations. 
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Neutron stars are not the sole objects in whose depth 
matter is compressed to high densities inaccessible in 
laboratory. Other representatives of the class of compact 
stars are white dwarfs and hypothetical quark stars [2]. 
While the size of a neutron star mainly depends on the 
balance between gravity force and degenerate neutron 
pressure, white dwarfs resist gravitational squeezing due 
to the electron degeneracy pressure, and quark or strange 
stars resist it due to the pressure of matter composed of 
quarks not combined into hadrons. Neutron stars are 
much more compact than white dwarfs. White dwarfs, 
regarded as a specific class of stars since the 1910s, with 
the mass M ~ Mq have the radius R ^ 10** km, which 
is comparable to Earth's radius but almost 1000 times 
greater than the radius of a neutron star [J]. There- 
fore, the matter density in their interiors is less than 
one-thousandth of pq. On the other hand, according to 
theoretical models, quark stars at M ~ Mq may be even 
more compact than neutron stars. But unlike neutron 
stars, quark stars have not been yet observed, and their 
very existence is questioned. 

Vitaly Lazarevich Ginzburg was among the pioneers 
of neutron star theoretical research. He predicted cer- 
tain important features of these objects before they were 
discovered by radio astronomers in 1967 and greatly con- 
tributed to the interpretation of observational data in the 
subsequent period. A few of his papers concerning these 
issues were published in 1964. In Refs. 4, 5] (the latter in 
co-authorship with L M Ozernoi), he described changes 
in the stellar magnetic field during collapse (catastrophic 
compression) and obtained the value B ~ 10^^ G, accu- 
rate to an order of magnitude for typical magnetic in- 
duction of a neutron star with a mass M ~ Mq , formed 
in the collapse. Moreover, Ginzburg derived expressions 
for the magnetic dipole field and the field uniform at 
infinity taking account of the space-time curvature near 
the collapsed star, in accordance with the general rela- 
tivity (GR). Today, these expressions are widely used 
to study magnetic neutron stars. In the same work, 
he predicted the existence of a neutron star magneto- 
sphere in which relativistic charged particles emit elec- 
tromagnetic waves in the radio to X-ray frequency range, 
and demonstrated the infiuence of magnetic pressure and 
magnetohydrodynamic instability and the possibility of 
detachment of the current-carrying envelope from the 
collapsing star [5]. Subsequent theoretical and observa- 
tional studies confirmed the importance of these prob- 
lems for the neutron star physics. In a one-and-a-half- 
page note \6] , Ginzburg and Kirzhnits formulated a num- 
ber of important propositions concerning neutron super- 
fiuidity in the interior of neutron stars (apparently inde- 
pendently of the earlier note by Migdal Q), the forma- 
tion of Feynman-Onsager vortices, a critical superfluid- 
ity temperature (Tc < 10^" K) and its dependence on 
the density {p ~ 10" - 10" g cm ), and the infiuence 
of neutron superfiuidity on heat capacity and therefore 
on the thermal evolution of a neutron star. These in- 
ferences, fully confirmed in later research, were further 
developed in Ginzburg's review [8*|, where he considered, 
inter alia, the superfluidity of neutrons and the supercon- 



ductivity of a proton admixture to the neutron fluid in 
the core of a neutron star. In Ref. [9], Ginzburg and Sy- 
rovatskii put forward the correct hypothesis of magnetic 
bremsstrahlung radiation from the source of X-rays in 
the Crab Nebula that turned out to be a plerion (a pul- 
sar wind nebula) surrounding the neutron star and of its 
origin from the envelope stripped off in collapse. In 1968, 
Ginzburg and co-workers greatly contributed to the elu- 
cidation of the nature of radio pulsars - cosmic sources 
of periodic radio pulses: they first developed the model 
of oscillating white dwarfs O, [ll[ , and later the mod- 
els of rotating neutron stars with strong magnetic fields 
[T2l - [T4j : in these studies, they discussed putative mech- 
anisms of pulsar radiation [isl . [T6| . Specific works were 
devoted to the role of pulsars in generation of cosmic rays 
[TtI I and the estimation of the work function needed to 
eject ions from the pulsar surface into the magnetosphere 
|18| . In 1971, Ginzburg published a comprehensive re- 
view [19], focused on the analysis of theoretical concepts 
of the neutron star physics and the physical nature of 
pulsars formulated by that time. The review contained 
a number of important original ideas, such as the estima- 
tion of typical neutron star magnetic fields B ^ 10^^ G, 
where he has pointed out that lower values to B ^ 10* G 
and higher values to B ~ 10^'^ - 10^^ G are also possible. 
These estimates were brilliantly confirmed by later stud- 
ies, which showed that a maximum in the magnetic field 
distribution of radio pulsars lies around B ^ 10^^ G 
[20l |: millisecond pulsars discovered in the 1980s have 
B ~ 10^-10" G [21] ; and fields of magnetars discov- 
ered in the 1990s reach up to B ~ lO^^ - 10" G [12. We 
must note that many of Ginzburg's other findings have 
wide application in neutron star research. Besides the fa- 
mous studies on superfiuidity and superconductivity, one 
should mention in this respect his investigations into the 
distribution of electromagnetic waves in the magnetized 
jjla sma, summarized in the comprehensive monograph 

Neutron stars are still insufficiently well known and 
remain puzzling objects despite their extensive studies 
in many research centers during the last 40 years. The 
aim of this review, from the perspective of modern astro- 
physics, is to highlight the main features of neutron stars 
making them unique cosmic bodies. It does not pretend 
to be comprehensive, bearing in mind the thousands of 
publications devoted to neutron stars, and is designed 
first and foremost to present the author's personal view 
of the problem. Fundamentals of neutron star physics 
are expounded in the excellent textbook by Shapiro and 
Teukolsky [24]. More detailed descriptions of develop- 
ments in selected branches of neutron star astrophysics 
can be found in monographs (such as [1,111]) and spe- 
cialized reviews (published, among others, in Physics- 
Uspekhi - i.e., [M, HI, 113] ) ■ 

In Sect. [2] we give general information on the physical 
properties of neutron stars and related physical and as- 
trophysical problems; the history of relevant research is 
also outlined. Section [3] equally concisely illustrates the 
"many faces" of neutron stars as viewed by a terrestrial 
observer. The following sections are of a less general 
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character, each being concerned with a specific aspect 
of astrophysics. The hst of these aspects is far from ex- 
haustive. For example, we will not consider the physics 
of the pulsar magnetosphere and the mechanisms under- 
lying generation of its radiation, which are dealt with in 
the voluminous literature (see, e.g., [H, [27l - l29t ). The 
same refers to nucleon superfluidity discussed in com- 
prehensive reviews (e.g., 26]). Only in passing are men- 
tioned the neutrino emission mechanisms, exhaustively 
described in the reviews of D G Yakovlev with coworkers 
[26l l30l | . The list of references on the problems concerned 
in this review is neither exhaustive, nor even representa- 
tive, which could not be otherwise, keeping in view the 
format of the article. The author apologizes to those re- 
searchers whose important contributions to the physics 
of neutron stars are not cited in this publication. A more 
complete (although still nonexhaustive) bibliography can 
be found in monograph [25j . 



2 Basic facts about neutron stars 



2.1 Neutron stars as relativistic objects 

Of great importance for neutron stars, unlike ordinary 
stars, are the GR effects [lH. The structure of nonro- 
tating stars is described by the nonrelativistic equation 
of hydrostatic equilibrium for a spherically symmetric 
body in GR, that is the Tolman-Oppenheimer-Volkojf 
(TOV) equation [S^, HI] . It also gives a very good ap- 
proximation for rotating stars, except those with mil- 
lisecond rotation periods. The minimal possible period 
is ~ 0.5 ms, but that observed to date is almost thrice 
as large, 1.396 ms 34], characteristic of the "slow ro- 
tation regime" in which the effects of rotation can be 
taken into account in terms of the perturbation theory 
[25I Ch. 6]. The corrections introduced by the magnetic 
field are negligibly small for the large-scale structure of 
a neutron star (at least for B < 10^^ G). The effects of 
the known magnetic fields B < 10^^ G can be important 
in stellar envelopes, as we will discuss in Sect. El Solu- 
tion of the TOV equation for a given equation of state of 
neutron star matter yields a family of stellar structure 
models, whose parameter is pc, the density in the center 
of the star. The stability condition requiring that M{pc) 
be an increasing function is satisfied within a certain 
range of stellar masses and radii; the maximum mass 
■M^max, compatible with the modern theory, is approxi- 
mately Minax ~ 1.5-2.5 Mq, depending on the equation 
of state being used, while the minimal possible mass of 
a neutron star is Mmin ~ 0.1 Mq. The significance of 
the GR effects for a concrete star is determined by the 
compactness parameter 

Xg = Tg/R, where Vg = 2GM/c^ sa 2.95 M/Mq km 

is the Schwarzschild radius, G is the gravitational con- 
stant, and c is the speed of light. Gravity at the stellar 



surface is determined by the equality 



9 = 



GM _ 1.328 X 10^4 M/Me 



® cm s-^ (2) 



where Rs = R/ilQ^ cm). 

The canonical neutron star is traditionally a star with 
M = 1.4 Mq and i? = 10 km (g = 2.425 x lO^^^ cm s^^). 
Note that the best and most detailed equations of state 
available to date predict a slightly lower compactness: 
i? « 12 km at M = 1.4 Mq (see Sect. g^). Substi- 
tuting these estimates in ([1]), we see that the effects of 
general relativity for a typical neutron star amount to 
tens of percent. This has two important consequences: 
first, the quantitative theory of neutron stars must be 
wholly relativistic; second, observations of neutron stars 
open up unique opportunities for measuring the effects 
of general relativity and verification of their prediction. 

The near-surface photon frequency (denoted by wq) in 
a locally incrtial reference frame undergoes a redshift to 
ujoo according to 
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Therefore, the thermal radiation spectrum of a star 
with an effective temperature Teg, measured by a dis- 
tant observer, is displaced toward longer wavelengths 
and corresponds to a lower effective temperature T^ = 
Tcs \/l - Xg. 

Along with the radius R, determined by the equato- 
rial length 2tiR in the locally incrtial reference frame, 
one often introduces the apparent radius for a distant 
observer: i?oo = ^ (1 + -s^g)- In particular, for the canon- 
ical neutron star we have Roo = 13 km, and for a more 
realistic model of the same mass, i?oo ~ 15 km. The ra- 
dius i? of a neutron star decreases as its mass increases, 
but the growth of Zg with the reduction of the radius 
and the increase in the mass leads to the appearance 
of a minimum in the dependence Roo{M). One can 
show that the apparent stellar radius cannot be smaller 
than i?™" = 7.66 (M/Mq) km (H]. The overall ap- 
parent photon luminosity oc Rl^{T^)'^ is related 
to the luminosity in the local stellar reference frame as 
L"^ = {1 ~ Xg) L^. The expressions for i?oo and are 
in excellent agreement with the notion of light bending 
and time dilation in the vicinity of a massive body. The 
light bending enables a distant observer to "look behind" 
the horizon of a neutron star. For example, the observer 
can simultaneously see both polar caps of a star having 
a dipole magnetic field at a proper dipole inclination an- 
gle to the line of sight. This effect actually occurs in 
observation of pulsars. Naturally, such effects must be 
taken into account when comparing theoretical models 
and observations. 

According to GR, a rotating star having a shape other 
than the ellipsoid of revolution can emit gravitational 
waves. Shape distortions may be caused by star oscil- 
lations and other factors. It has been speculated [111 
that gravitational waves emitted by rapidly rotating neu- 
tron stars can be recorded by modern gravitational an- 
tennas. However, these antennas appear more suitable 
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for recor ding gravitational waves from merging neutron 
stars (36l [37|. 

While gravitational waves have not yet been detected 
by ground-based antennas, they have already been regis- 
tered in observations of "space antennas" - double neu- 
tron stars. Two stars orbiting a common center of masses 
are known to emit gravitational waves. The first pulsar 
rotating in an orbit together with another neutron star 
was discovered by Hulse and Taylor in 1974 (Nobel Prize 
of 1993). It remained the sole known object of this kind 
for over 20 years. At least 9 such systems have been 
described to date. The most remarkable of them is the 
double pulsar system J0737-3039, a binary in which both 
neutron stars are seen as radio pulsars [2ll . ISSj . 

The known binary systems of neutron stars have com- 
pact orbits and short periods of revolution. The orbital 
period of the Hulse -Taylor pulsar is less than 8 hours, 
and the large semiaxis of the orbit is about two million 
kilometers, or almost two orders of magnitude smaller 
than the distance between the Sun and Earth. Gravi- 
tational radiation is so strong that the loss of energy it 
carries away results in a significant decrease in both the 
orbit size and the orbital period. The measured decrease 
in the orbital period of the Hulse - Taylor pulsar is con- 
sistent with that predicted by GR within a measurement 
error to a few tenths of a percent. 

Another GR effect is the periastron shift or relativistic 
precession of the orbit that is orders of magnitude greater 
than Mercury's perihelion shift (to be precise, those 7.5% 
of its perihelion shift, namely 0.43" per year, that can- 
not be accounted for by the influence of other objects of 
the Solar system are explained in the GR framework). 
For example, for the Hulse -Taylor pulsar the relativis- 
tic periastron shift is 4.22° per year, and for the double 
pulsar it is 16.9° per year. 

The third measured effect is geodesic precession of a 
rotating body that moves in an orbit, a precession anal- 
ogous to the spin-orbital interaction in atomic physics. 
The measurement of geodesic precession made it possi- 
ble to reconstruct the time dependence of the direction 
of the Hulse -Taylor pulsar magnetic axis. It turned out 
that its directivity pattern would no longer intersect the 
line of sight of an Earth-based observer around 2025, 
and the pulsar would become invisible for two centuries 
[39, 40]. One cannot exclude that the companion star 
will become visible. 

The double pulsar proved an even better laboratory 
than the Hulse - Taylor pulsar for the verification of the 
effects of general relativity. First, registration of radio 
pulses from both neutron stars of the binary system al- 
lows directly measuring the radial velocities from the 
Doppler shift and geodesic precession of either star from 
the altered pulse shape. Second, the line of sight of a ter- 
restrial observer lies virtually in the plane of the double 
pulsar orbit (at the inclination to the normal « 89°). 
This permitted for the first time to reliably measure 
the so-called Shapiro delay parameters (two parameters 
characterizing the time delay of an electromagnetic wave 




Figure 1: Schematic of the neutron star structure, from 
the paper [11] by V L Ginzburg. The labels read (from 
the center outwards): Core, Superfiuid neutron liq- 
uid. Superconducting proton liquid. Rigid crust. Dense 
plasma envelope. 

passing a starQ Five of the seven independent post- 
Kepler parameters characterizing the effects of general 
relativity were measured for the double pulsar. Any two 
of them uniquely define masses of both pulsars Ma and 
Mb , while the measurement of the remaining ones may 
be regarded as verification of GR. HR has brilliantly 
passed this test: any two of the measured parameters 
gave the same values Ma = 1.337 Mq, Mr = 1.249 Mg 
within measurement errors < 0.001 |38| . 

One more unique relativistic object discovered in 2005 
is the pulsar PSR J1903-I-0327, having the rotation pe- 
riod 2.15 ms and moving in an inclined highly elliptical 
orbit (eccentricity e = 0.44 and inclination 78°) in a 
pair with a main-sequence star (an ordinary star of a 
mass M w Mq, identified in the infrared [41]). Observa- 
tions with the Arecibo radio telescope yielded three post- 
Kepler parameters: both Shapiro delay parameters and 
orbital precession. The estimate M — 1.67 ± 0.01 Afg 
was obtained for the pulsar mass under the assumption 
that the orbital precession is due to the effects of GR 
alone 02] ■ It is the largest mass of neutron stars mea- 
sured thus far. We note, however, that the influence 
of nonrelativistic effects, such as tides at the compan- 
ion star caused by gravitational attraction of the pul- 
sar, on the orbital precession cannot be totally excluded. 
With this uncertainty, a more conservative estimate is 
M = 1.67±0.11 Mq [22]. (See the Note added in proof). 

2.2 The biggest enigmas of the neutron 
star structure 

Two main qualitatively different regions, the core and 
the envelope, are distinguished in a neutron star. The 
core is in turn subdivided into the outer and inner core, 
and the envelope into the solid crust and the liquid 

^In addition, the proximity of the orbit plane to the line of 
sight allowed one to observe modulation of pulsed emission from 
one pulsar passing through the atmosphere of the other, which 
provided supplementary information on their magnetic fields and 
magnetospheres . 
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ocean. Such a division into four essentially different re- 
gions was proposed in the review by Ginzburg of 1971 
UM (Fig- [U- The outer core of a neutron star is usu- 
ally several kilometers thick, and its matter density is 

0. 5po ^ P ^ 2po- This matter, accounting for the largest 
fraction of the stellar mass, has well-known qualitative 
characteristics (see, e.g., [H, Ch. 5 and 6]). It is a neu- 
tron superfluid with an admixture of the superconduct- 
ing proton component (see Fig. [T]), as well as electrons 
and /z~-mesons (muons), all these constituents being 
strongly degenerate. The inner core occupies the central 
part with p ^ 2po and has a radius to several kilome- 
ters. It can be present in rather massive neutron stars, 
M > 1.4- 1.5 M0 (in less massive neutron stars, density 
does not reach 2po)- Neither the composition nor the 
properties of matter in the inner core are known because 
the results of their calculation strongly depend on the 
theoretical description of collective fundamental interac- 
tions. From this standpoint, studies of neutron stars are 
important not only for astrophysics but also for nuclear 
and elementary particle physics. The available theoreti- 
cal models presume the following hypothetical options: 

1. hyperonization of matter - the appearance of vari- 
ous hyperons (first of all. A- and S^-hyperons); 

2. pion condensation - formation of a Bose condensate 
from collective interactions with the properties of 
TT-mesons; 

3. kaon condensation - formation of a similar conden- 
sate from if -mesons; 

4. deconfinement - phase transition to quark matter. 

The last three variants, unlike the first one, are not 
feasible for all modern theoretical models of matter 
with supranuclear density; therefore, they are frequently 
called exotic [H Ch. 7]. 

According to current concepts, the core of a neutron 
star contains superfluid baryonic matter. Superfluidity 
reduces the heat capacity of this matter and the neutrino 
reaction rate. However, superfluidity may be responsi- 
ble for an additional neutrino emission due to Cooper 
pairing of nucleons at a certain cooling stage in those 
parts of the star where the temperature decreases below 
critical values. These effects and their influence on the 
cooling rate of a neutron star are reviewed in [2^ . 

In the stellar envelopes the matter is not so extraor- 
dinary: the atomic nuclei are present separately there. 
Nevertheless, this matter also occurs under extreme 
(from the standpoint of terrestrial physics) conditions 
that cannot be reproduced in the laboratory. This makes 
such matter a very interesting subject of plasma physics 
research P]. Equally important is the fact that an ade- 
quate theoretical description of stellar envelopes is indis- 
pensable for the correct interpretation of characteristics 
of the electromagnetic radiation coming from the star, 

1. e., for the study of its core by means of comparison of 
theoretical models and astronomical observations. 



2.3 The birth, hfe, and death of a neu- 
tron star 

A neutron star is a possible end product of a main- 
sequence star ("normal" star) Neutron stars are 
believed to be formed in type-I supernovae explosions 
[43| - |47| . An explosion occurs after a precursor to a su- 
pernova has burned out its nuclear "fuel" : flrst hydro- 
gen, then helium produced from hydrogen, and finally 
heavier chemical elements, including oxygen and mag- 
nesium. The end product of subsequent nuclear trans- 
formations is isotopes of iron-group elements accumu- 
lated in the center of the star. The pressure of the 
electron Fermi gas is the sole factor that prevents col- 
lapse of such an iron-nickel core to its center under the 
force of gravity. But as soon as a few days after oxygen 
burning, the mass of the iron core increases above the 
Chandrasekhar limit equal to 1.44 Mq, which is the max- 
imum mass whose gravitational compression is still coun- 
teracted by the pressure of degenerate electrons. Then 
gravitational collapse, i.e., a catastrophic breakdown of 
the stellar core, occurs. It is accompanied by the liber- 
ation of an enormous gravitational energy (> 10^^ erg) 
and a shock wave that strips off the outer envelopes of 
the giant star at a speed amounting to 10% of the speed 
of light, while the inner part of the star continues to 
contract at approximately the same rate. The atomic 
nuclei fuse into a single giant nucleus. If its mass sur- 
passes the Oppenheimer - Volkoff limit, that is the maxi- 
mum mass that the pressure of degenerate neutrons and 
other hadrons is able to support against gravitational 
compression (« 2-3 M©), according to modern theoreti- 
cal models) , the compression cannot be stopped and the 
star collapses to form a black hole. It is believed that 
the collapse resulting in a black hole may be responsi- 
ble for the flare from a hypernova, hundreds of times 
brighter than a supernova; it may be a source of mys- 
terious gamma-ray bursts coming from remote galaxies 
[48*, "4^ . If the mass remains below the Oppenheimer - 
Volkoff limit,, a neutron star is born whose gravitational 
squeezing is prevented by the pressure of nuclear matter. 
In this case, about 1% of the released energy transforms 
into the kinetic energy of the envelopes flying apart, 
which later give rise to a nebula (supernova remnant), 
and only 0.01% (~ 10^® erg) into electromagnetic radi- 
ation, which nevertheless may overshine the luminosity 
of the entire galaxy and is seen as a supernova. 

Not every star completes its evolution as a supernova 
(not to mention a hypernova); only massive stars with 
M > 8 M0are destined to have such a fate. A less mas- 
sive star at the end of its lifetime goes through a giant 
phase, gradually throwing off the outer envelopes, and 
its central part shrinks into a white dwarf. 

A newborn neutron star has the temperature above 
10^° -10^^ K ; thereafter, it cools down (rather fast ini- 
tially, but slower and slower afterwards), releasing the 
energy in the form of neutrino emission from its depth 
and electromagnetic radiation from the surface. But the 
evolution of a neutron star is not reduced to cooling 
alone. Many neutron stars have strong magnetic fields 
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that also evolve through changes in strength and config- 
uration. A rotating neutron star having a strong mag- 
netic field is surrounded by an extended plasma magne- 
tosphere formed due to the knockout of charged particles 
from the surface by the rotation-induced electric field, 
thermal emission, and the birth of electron-positron pairs 
upon collisions of charged particles of the magnetosphere 
with one another and with photons. Given a sufficiently 
fast rotation of a star, its magnetosphere undergoes col- 
lective acceleration of the constituent particles in the 
parts where plasma density is too low to screen the 
strong electric field induced by rotation. Such processes 
generate coherent directed radio-frequency emission due 
to which the neutron star can be seen as a radio pulsar if 
it rotates such that its directivity pattern intersects ob- 
server's line of sight. The rotational energy is gradually 
depleted and the particles born in the magnetosphere 
have a charge whose sign is such that the induced elec- 
tric field makes them propagate toward the star; they 
accelerate along the magnetic force lines, hit the star 
surface near its magnetic poles, and heat these regions. 
A similar process of heating magnetic poles occurs in 
the case of accretion (infall of matter) onto a star, e.g., 
as it passes through dense interstellar clouds or as the 
matter outfiows from the companion star in a binary sys- 
tem. The hot polar caps emit much more intense X-rays 
than the remaining surface; as a result, such neutron 
stars look like X-ray pulsars. Pulsed X-ray radiation is 
also observed from thermonuclear explosions of accreted 
matter at the surface of a rotating neutron star (see, e.g., 
review [50l|). 

Cooling, changes in the magnetic field, or a slowdown 
of rotation may cause starquakes, associated with varia- 
tions of the crustal shape, phase transformations in the 
core, and interaction between the normal and superfiuid 
components of the core and the crust [5ll - l56l | . Starquakes 
are accompanied by liberation of the thermal energy and 
sharp changes in the character of rotation [S^]. More- 
over, the matter falling onto the star during accretion 
undergoes nuclear transformations at the surface and, 
due to its weight, causes additional transformations in 
the depth of the envelope that alter the nuclear com- 
position and liberate energy [SSl-GO] . In other words, 
neutron stars not only cool down but also are heated 
from the inside. 

A single neutron star eventually exhausts its supply of 
thermal and magnetic energy and fades away. A star has 
more promising prospects for the future if it is a member 
of a binary system. For example, if the companion star 
overfills its Roche lobe (the region in which matter is 
gravitationally coupled to the companion), this matter 
accretes onto the neutron star so intensely as to make it a 
bright source of X-ray radiation by virtue of the released 
gravitational and thermonuclear energy. In this case, 
the infiowing matter forms an accretion disk around the 
neutron star, which also radiates X-ray emission, and 
this luminosity changes with time, e.g., as a result of 
disk precession or variations of the accretion rate. The 
character of accretion strongly depends on neutron star 
magnetization and the rotation period [61.] . If the mass 



of the accreted matter surpasses a critical threshold, the 
neutron star collapses into a black hole. 

If the companion of a neutron star is a compact ob- 
ject, then the radius of their mutual orbit may be small 
enough to enable gravitational waves emitted by such a 
system to appreciably infiuence its evolution. The or- 
bital radius of the compact binary system decreases as 
gravitational radiation continues until the two compan- 
ions merge together, giving rise to a black hole with the 
release of enormous gravitational energy comparable to 
the stellar rest energy Mc^ ~ 10^^ erg in the form of 
neutrino and gravitational radiation (this will happen to 
the Hulse - Taylor pulsar and the double pulsar in some 
300 and 85 million years, respectively). 

One can say that the three main driving forces of the 
evolution of a neutron star, responsible for its observa- 
tional manifestations, are rotation, accretion, and mag- 
netic field. 

2.4 The formation of neutron star con- 
cepts 

Baade and Zwicky [G^ theoretically predicted neutron 
stars as a probable result of supernova explosions less 
than 2 years after the discovery of the neutron [63l|.Thev 
also put forward the hypothesis (now universally ac- 
cepted) that supernovae are important sources of galac- 
tic cosmic rays and coined the term "supernova" itself 
to differentiate between these unusually bright objects 
formed in a gravitational collapse giving rise to a neu- 
tron star from more numerous nova stars originating, 
as known today, from thermonuclear burning of the ac- 
creted matter at the surface of white dwarfs. The pop- 
ular belief that Landau predicted neutron stars in 1932 
[131, based on recollections by Leon Rosenfeld |6J], is 
not accurate:the meeting of Landau with N Bohr and 
L Rosenfeld occurred in 1931, before the discovery of 
neutrons. Nevertheless, it is true that Landau already 
foresaw the existence of neutron stars at that time and 
suggested a hypothesis according to which stars with a 
mass greater than 1.5 Mq have a region in their inte- 
rior where the density of matter "becomes so great that 
atomic nuclei come in close contact, forming one gigantic 
nucleus" p5| . 

Not a single neutron star was observed for 43 years 
after Baade and Zwicky's prediction. But theorists con- 
tinued to work. In 1938, Zwicky [6^ estimated the max- 
imum binding energy of a neutron star and the gravi- 
tational red shift of photons emitted from its surface. 
A few months later, Tolman [32] and Oppenheimer and 
Volkoff [33| derived the aforementioned TOV equation; 
moreover, the latter authors computed the limiting mass 
of a neutron star, Mmax, although it proved underesti- 
mated because they ignored baryon-baryon interactions. 
Equations of state of nuclear matter began to be exten- 
sively studied in the 1950s. In 1959, Cameron [63] ob- 
tained the first realistic estimate of Mmax ~ 2 Mq . He 
was the first to show that the core of a neutron star may 
contain hyperons. In the same year, Migdal fzl], based 
on the concept of superfluidity in atomic nuclei proposed 
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by A Bohr, Mottelson, and Pines [6l|, predicted super- 
fluidity of neutron star matter. In 1960, Ambartsumyan 
and Saakyan p9| constructed tlie equation of state of 
superdense matter by taking electrons, muons, and all 
hadrons known at that time into consideration. They 
hypothesized that the core of a neutron star consists of 
two components, the outer composed of nucleons and 
the inner containing hyperons. The following year, Zel- 
dovich [T^I derived an extremely stiff equation of state 
of a neutron star in which the speed of sound tends to 
the speed of light as the density increases. Finally, in 
the 1960s, the first estimates of neutrino emission from 
the interior of a neutron star J2j and its cooling 
[76j | were reported; in additiori, the presence of a strong 
magnetic field was predicted [j, Q and the deceleration 
of rotation of a magnetized neutron star due to magne- 
todipole radiation calculated fn\ . 

The first simplest models of neutron star cooling al- 
ready demonstrated that the surface temperature of a 
typical neutron star might be as high as hundreds of 
thousands or millions of degrees, meaning that the star 
emits thermal radiation largely in the form of soft X-rays 
unable to penetrate Earth's atmosphere. The progress 
in astronautics and the advent of X-ray astronomy in 
the early 1960s ^] gave hope that such radiation would 
be detected in outer space. However, it took almost 30 
years to reliably identify thermal X-ray components in 
the spectra of neutron stars by the X-ray telescope on 
board the ROSAT satellite, witch produced images with 
a resolution of a few angular seconds [79|] . 

Other means were proposed to search for neutron 
stars. Zeldovich and Guseinov [S^ suggested that they 
could be detected in binary systems with optical stars 
from the Doppler shifts of optical spectral lines. Karda- 
shev [U and Pacini [77J put forward the correct hy- 
pothesis that the rotational energy of a neutron star 
was transferred via the magnetic field to the surround- 
ing nebula formed in the collapse at stellar birth. These 
authors regarded the Crab Nebula as such a candidate 
object. However, neutron stars unexpectedly manifested 
themselves as radio pulsars. 

As known, the first pulsar was discovered by radio 
astronomers at Cambridge in 1967 [82| (a retrospective 
study of the archives of the Cambridge group gave evi- 
dence of radio pulses from pulsars dating back to 1962- 
1965 ^3]). Antony Hewish, who headed the research 
team, was awarded the Nobel Prize in physics in 1974 
for this achievement. A correct explanation of these ob- 
servations soon after their publication was proposed by 
Thomas Gold in the paper entitled "Rotating neutron 
stars as the origin of the pulsating radio sources" [84| . 
It is less widely known that Shklovsky [SSj arrived at 
the conclusion, based on analysis of X-ray and optical 
observations, that radiation from Scorpio X-1 (the first 
X-ray source discovered outside the Solar System 78])) 
originated from the accretion of matter onto a neutron 
star from its companion. Unfortunately, this conclusion 
(later fully confirmed [S^) was accepted too sceptically 
at that time [87(|. 

The discovery of pulsars gave powerful impetus to the 



development of theoretical and observational studies of 
neutron stars. With over one thousand publications de- 
voted to these celestial bodies appearing annually, a new 
class of astronomical objects containing neutron stars is 
discovered once every few years. For example. X-ray 
pulsars were described in 1971, bursters (sources of X- 
ray bursts) in 1975, soft gamma repeaters (SGR5) in 
1979, millisecond pulsars in 1982, radio-silent neutron 
stars in 1996, anomalous X-ray pulsars (AXPs) in 1998, 
and rapid (rotating) radio transients (RRATs) in 2006. 
Clearly, it is impossible to cover all these developments 
in a single review. We try instead to depict the cur- 
rent situation in certain important branches of the the- 
ory, although we start with observable manifestations of 
neutron stars. 

3 Observational manifestations 
of neutron stars 

Radiation from neutron stars is observed in all ranges of 
the electromagnetic spectrum. As well as 40 years ago, 
most of them (about 1900 as of 2010 [13]) are seen as 
radio pulsars. Some 150 of the known neutron stars are 
members of binary systems with accretion and manifest 
themselves largely in the form of X-ray radiation from 
the accretion disk or flares produced by explosive ther- 
monuclear burning in the star outer layers. Certain such 
systems make up X-ray transients in which periods of ac- 
tive accretion (usually as long as several days or weeks) 
alternate with longer periods of quiescence (months or 
sometimes years) during which X-ray radiation from the 
hot star surface is recorded. In addition, over one hun- 
dred isolated neutron stars are known to emit X-ray ra- 
diation. 

3.1 Cooling neutron stars 

A large fraction of emission from isolated neutron stars 
and X-ray transients in quiescence appears to originate 
at their surface. To interpret this radiation, it is very 
important to know the properties of envelopes contribut- 
ing to the spectrum formation. Conversely, comparison 
of predictions and observations may be used to deduce 
these properties and to verify theoretical models of dense 
magnetized plasma. Moreover, investigating the proper- 
ties of the envelopes provide knowledge of the parameters 
of a star as a whole and of the observational constraints 
on such models. 

For each theoretical model of neutron stars, a cooling 
curve describes the dependence of the overall photon lu- 
minosity in the frame of a distant observer on time 
t elapsed after the birth of the star (see review [Hi and 
the references therein). 

There are not many neutron stars in whose spectrum 
the cooling-related thermal component can be distin- 
guished from the emission produced by processes other 
than surface heating, e.g., those proceeding in the pul- 
sar magnetosphere, pulsar nebula, and accretion disk. 
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Fortunately, there are exceptions [89i , such as relatively 
young {t < lO^ears) pulsars J1119-6127, B1706-44, 
and Vela, whose spectra are readily divisible into ther- 
mal and nonthermal components, and medium-aged (t ^ 
10<^ years) pulsars B0656-hl4, B1055-52, and Geminga. 
The spectra of the latter three objects, dubbed "three 
musketeers" [zl, are fairly well described by the 
three-component model (power-law spectrum of mag- 
netospheric origin, thermal spectrum of the hot polar 
caps, and thermal spectrum of the remaining surface). 
Even more imp ortant is the discovery of radio-silent neu- 
tron stars [90j,with purely thermal spectra. These are 
central compact objects (CCOs) in supernova remnants 
91 1 and X-ray "dim" isolated neutron stars (XDINS) 
9^ . Observations indicate that CCO may have mag- 
netic fields B ~ 10^° - 10" G (shghtly weaker than in the 
majority of normal pulsars but stronger than in millisec- 
ond pulsars) , and XDINSs may have magnetic fields 
B > 10^'^ G (somewhat stronger than ordinary) [23] ■ As 
many as 7 XDINSs are constantly known during the last 
decade, and astrophysicists call them the "Magnificent 
Seven" [H, IH, [Hi . By the way, confirmed CCOs count 
also seven, but three more objects are candidates wait- 
ing to be included in the list [93] . The spectra of at least 
5 radio-quiet neutron stars exhibit wide absorption lines 
for which no fully satisfactory theoretical explanation 
has been proposed thus far. Certain authors hypothe- 
size that they can be attributed to ionic cyclotron har- 
monics in a strong magnetic field, but rigorous quantum 
mechanical calculations |95i , i96i] have proved that such 
harmonics in neutron star atmospheres are too weak to 
be observed. 

Besides cooling processes, heating processes of differ- 
ent natures occur in neutron stars. Sometimes they com- 
pete with the heat delivered to the stellar surface from 
the core and must be taken into consideration. Such 
stars include: 

— old neutron stars {t > 10® years) for which the cool- 
ing curves would go down to the low temperature 
region (Toff < 10^ K) when disregarding heating; 

— magnetars, which are relatively young (t < 10* 
years) neutron stars with superstrong {B > 10^^ G) 
magnetic fields manifested as AXPs and SGRs 
[13, [131 • The strong X-ray luminosity of magne- 
tars cannot be explained by the "standard cooling 
curve" [11] . Thompson [97] suggested that it should 
be ascribed to heating due to dissipation of a su- 
perstrong magnetic field. Recent studies (98l - ll00l ] 
provided some support to this hypothesis. 

The discussion of neutron star cooling is to be con- 
tinued in Section [T] For now, we emphasize that the 
processes of cooling, heating, and heat transfer turn the 
surface of a neutron star into a source of thermal radia- 
tion with a spectral maximum in the soft X-ray region. 

3.2 Pulsars 

Pulsed radiation related to the proper rotation of neu- 
tron stars contains important additional information 



[2^. [28l. [lOlj . For example, simultaneous measurement 
of emission at several radio frequencies allows determin- 
ing the dispersion measure from the phase shift, which in 
turn permits roughly estimating the distance to a pulsar. 
Measurements of the pulsation amplitude of the thermal 
component in the spectrum characterize the nonuniform 
temperature distribution over the surface. The period 
of pulsations P and its time derivative P for isolated 
(nonaccreting) pulsars give an idea of the star magnetic 
field (of its dipole constituent, to be precise) and age: 

lO^^-^y^PP/l s G, t - ipsR EE 0.5P/P, (4) 

where tpsR is the so called characteristic age of a pul- 
sar. Such estimate makes no sense for accreting pulsars 
because their period of rotation may depend on the in- 
teraction between the magnetic field and the accretion 
disk. For example, the rotation can be accelerated by 
virtue of the transfer of angular momentum from matter 
falling onto the pulsar; then P < 0. 

We note that the age of a neutron star can be deduced 
from the age of the remnant of the supernova hosting this 
star. As a rule, the age of the remnant is inferred with 
an error 10% from the rate at which the envelopes 
fly apart from each other. Certainly, the remnant it- 
self must be accessible to observation, which is rarely 
the case (as a rule, only a.t t < 10'* years). When both 
the characteristic age and the age of the remnant are 
known, they are consistent with each other to the order 
of magnitude, but their numerical difference may be as 
great as 2-3-fold. This means that none of the age esti- 
mation methods is entirely reliable. The exceptions are 
5 superno vae w hose flares are documented in historical 
chronicles |102| . 

Simultaneous measurement of the age and magnetic 
field permits imposing limitations on the decay rate of 
the stellar magnetic field. The relevant theoretical esti- 
mates are significantly different, depending on the field 
configuration inside the crust and the core, which in turn 
depends on the model and on the hypothesis of the na- 
ture of the field; they also depend on the electric con- 
ductivity included in a given model, which is in turn 
a function of the poorly known chemical composition of 
the envelopes, the microscopic structure of the crust, the 
content of admixtures, and the defects of the crystal lat- 
tice. 

X-ray radiation from pulsars, similarly to radio emis- 
sion, carry important information. Generally speaking, 
the X-ray spectrum of pulsars contains both thermal 
and nonthermal components. The latter is produced in 
the magnetosphere either by synchrotron radiation or 
inverse Compton scattering of charged particles acceler- 
ated to relativistic energies by magnetospheric electro- 
magnetic fields. This component is usually described by 
a power-law spectrum. The thermal component is di- 
vided into "hard" and "soft" constituents. The former 
is supposed to be a result of radiation from the polar 
caps heated to millions of degrees, where the magnetic 
field does not substantially deviate from the normal to 
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the surface. In the dipole field model, the radius of these 
regions is estimated as 

1 /2 

i?cap « (^) « 0.145i?^/^ (P/1 s)-V2 km. (5) 

The soft constituent corresponds to radiation from the 
remaining, cooler surface that may be due to the heat 
coming up from the stellar core or inner crust. 

3.2.1 Ordinary pulsars 

Ordinary pulsars are isolated pulsars with periods from 
tens of milliseconds to several seconds. Their character- 
istic magnetic field given by formula ([4]) varies from a 
few gigagauss to 10"^'* G with typical values B ^ 10^^ — 
10^'^ G, and the characteristic age is from several cen- 
turies to 10^" years with typical values ipsR ^ 10^ — 10^ 
years [20l |. The X-ray spectrum of thermal radiation 
from certain normal pulsars has been measured, which 
allows the cooling theory methods to be applied to their 
study. Unlike the case of neutron stars lacking pulsation, 
the estimates of ipsR and B make the class of thermal 
structure models more definite and thereby restrict the 
scattering of possible cooling curves. 

3.2.2 Millisecond pulsars 

Millisecond pulsars have magnetic fields B ^ 10^ — 
10^*^ G and are from tens of millions to hundreds of bil- 
lions of years old, with typical values ipsR ^ 10^ — 10^*^ 
years (with the exception of PSR J0537— 6910 with ab- 
normally high P) [23|. The relatively weak magnetic 
field and the short period of millisecond pulsars may re- 
sult from the pulsar having passed through the stage 
of accretion in the course of its evolution, which re- 
duced the magnetic field and increased the angular mo- 
mentum due to the interaction between the accreting 
matter and the magnetic field [2ll, [6l|. Certain iso- 
lated millisecond pulsars emitting in the X-ray range 
show a thermal constituent in their s pect ra produced 
by radiation from the hot polar caps |103j . It is con- 
venient to write formula ([5]) for millisecond pulsars in 

the form i?cap ~ Rq^^ (-F'/21 ms)~^/^ km, which readily 
shows that the hot region covers a large surface area of 
the pulsar. 

3.2.3 Anomalous pulsars 

Many neutron stars manifest themselves through pulsed 
radiation in the X-ray part of the spectrum. They are 
referred to as X-ray pulsars. Some of them are located 
in binary stellar systems. Evidently, those radio pulsars 
that exhibit thermal radiation from polar caps are also 
X-ray pulsars. 

Unlike these "normal" X-ray pulsars, AXPs have an 
unusually long period, P sa 6 — 12 s and high X-ray lu- 
minosity 10^3 — iq35 gj.g g-i^ being in the same time 
isolated [2^ . |104| . Their magnetic fields and character- 
istic ages estimated from (j4]), suggest that these objects 



may be magnetars. An alternative explanation of their 
properties is based on the assumption that they are neu- 
tron stars with "normal" magnetic fields B ~ 10^^ G 
that slowly accrete matter fr om the disk remaining af- 
ter the supernova explosion [105l |. In other words, the 
nature of these objects remains obscure. 

3.3 Neutron stars in binary systems 

A neutron star in a binary system is paired with another 
neutron star, a white dwarf, or an ordinary (nondegen- 
erate) star. Binary systems containing a neutron star 
and a companion black hole are unknown. Measuring 
the parameters of the binary orbit provides additional 
characteristics of the neutron star, e.g., its mass. 

The infall of matter onto a neutron star is accompa- 
nied by the liberation of energy, which turns the system 
into a source of bright X-ray radiation. Such systems are 
categorized into markedly different subclasses: low-mass 
X-ray binary systems with a dwarf (either white or red) 
of mass < 2Mq as the companion, and relatively short- 
lived massive systems in which the mass of the compan- 
ion star is several or tens of times greater than AfQand 
accretion of matter onto the neutron star is extremely 
intense. 

X-ray binaries may be sources of regular (periodic) and 
irregular radiation, and are subdivided into permanent 
and temporary (transient). Emission from some of them 
is modulated by neutron star rotation, others are sources 
of quasiperiodic oscillations (QPOs), bursters (neutron 
stars whose surface from time to time undergoes explo- 
sive thermonuclear burning of the accreted matter) , and 
so on. 

The QPOs, first observed in 1985 [106j . occur in X-ray 
binary systems containing compact objects, such as neu- 
tron stars (typically within low-massive X-ray binaries), 
white dwarfs, and black holes. There are a variety of 
hypotheses on the nature of the QPOs (see 107 109]). 
They seem to originate in the accretion disk. According 
to some hypotheses, they are related to the Kepler fre- 
quency of the innermost stable orbit permitted by GR, a 
resonance in the disk, or a combination of these frequen- 
cies with the rotation frequency of the compact object. 
Given that one of these hypotheses is true, the QPOs in 
low-massive X-ray binary systems may become a tool for 
determining the parameters of neutron stars. 

X-ray luminosity of type-I bursts in bursters may 
reach the Eddington limit L^dd ~ 1-3 x 10^^ (M/Mq) 
erg s~^, at which radiation pressure on the plasma due 
to the Thomson scattering exceeds the force of gravity. 
Such flares are of special interest in that simulation of 
their spectra and intensity permits estimating the pa- 
rameters of a neutron star |110L lllll | . 

The spectra of certain soft X-ray transients during 
"quiescent" periods exhibit the thermal radiation com- 
ponent of the neutron star located in the system This 
allows comparing the cooling curves with observations, 
as in the case of isolated neutron stars, with the sole 
difference that the energy release due to accretion must 
be taken into account. On the one hand, this introduces 
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an uncertainty in the model, but on the other hand it 
permits verifying theoretical considerations concerning 
accretion onto the neutron star and thermonuclear trans- 
formations of matter in its envelopes. 

Of special interest are quasipermanent transients, i.e., 
those whose active and quiescent periods last a few y ears 
or longer. According to the model proposed in [112| . the 
thermal radiation in the periods of quiescence is due to 
the crust cooling after deep heating by accretion. Such 
cooling is independent of the details of the star structure 
and composition and therefore its analysis directly yields 
information on the physics of envelopes. Three sources 
of this kind are known: KS 1731-260, MXB 1659-29, 
and AX J1754.2-2754 [lll|. The crust of a neutron 
star gets heated during a long period of activity and re- 
laxes to the quasiequilibrium state in the quiescence pe- 
riod that follows, meaning that evolution of the thermal 
spectrum contains information about the properties of 
the crust. Therefore, analysis of the thermal luminosity 
of such an object and its time dependence provides in- 
formation about heat capacity and thermal conductivity 
of the crust in the period of activity preceding relaxation 
and about the equilibrium luminosity at rest. Such data 
may in tur n be used to obtain characteristics of the star 
as a whole |114j . 

Estimation of neutron star masses from the Kepler 
parameters of X-ray binary systems is not yet a reli- 
able method because of theoretical uncertainties, such as 
those related to the transfer of angular momentum via 
accretion. The most accurate estimates are obtained for 
binary systems of two neutron stars due not only to the 
absence of accretion but also to the marked GR effects, 
whose measurement allows determining the complete set 
of orbital parameters. Sufficiently accurate mass esti- 
mates (with an error < 0.2 M0) are also available for 
several binary systems containing a white dwarf as the 
companion and for the system mentioned in Section 12. 1[ 
in which the companion is a main-sequence star. These 
estimates lie in the range 1.1 Mq < M < 1.7 Mq. 

4 The core of a neutral star and 
supranuclear density matter 

In this section, we focus on the equation of state of the 
neutron star core, disregarding details of the microscopic 
theory and nonstationary processes. We note that the 
Fermi energy of all particles essential for this equation is 
many orders of magnitude higher than the kinetic ther- 
mal energy. Therefore, a good approximation is given by 
the equation of state of cold nuclear matter in which the 
dependence of the pressure on density and temperature, 
P{p,T) , is replaced by a one-parametric dependence 
P{p) at T ^ 0. 

4.1 The outer core 

Nucleons in the outer core of a neutron star form a 
strongly interacting Fermi liquid, whereas leptons make 



up an almost ideal Fermi gas. Therefore, the energy 
density £ can be represented as the sum of three terms, 

£{nn,np,ne,nf_,) = £Ar(n„,np) -I- fe("-e) (6) 

where rie, n^, n„, and Up are concentrations of electrons, 
/i~-mesons, neutrons, and protons. The equations of 
state and concentration of particles are determined by 
the energy density minimum at a fixed baryon volume 
density Ub ~ rin + rip and under the electroneutrality 
condition rie -I- = rip. This implies that the relations 
P'Ti — fJ'p + Me and /i^ = pe for chemical potentials /ij of 
particles j = ri,p, e,/i~, expressing the equilibrium con- 
ditions with respect to the electron and muon beta-decay 
and beta-capture reactions: n— >p + e + i7e, p + e ^ 
n + n p + p, + Vfj_, and p + fi ^ n + Ufj,, where 
i^cfi and Ve^p. are electron and muon neutrinos and an- 
tineutrinos. Neutron star matter (unlike the matter of a 
protoneutron star, i.e., the collapsed core within the first 
minutes after the supernova explosion) is transparent to 
neutrinos: therefore, the chemical potentials of neutrino 
and antineutrino are equal to zero. Electrons at the den- 
sities being considered are ultrarelativistic, and therefore 
Me ~ cpFe ~ 122.1 (rie/O.OSno)^/^ MeV, where ppe is the 
electron Fermi momentum, and uq — 0.16 fm~^ is the 
normal nuclear number density, which corresponds to 
the normal nuclear mass density po- In the general case, 
muons are moderately relativistic, which dictates the use 

of the general expression fj,^ = m^c^^l + p|^/(?7i^c)^. 

As soon as the equilibrium is known, the pressure can 
be found from the equation P = n^d(f /n{,)/dnb. 

Thus, the construction of the equation of state for the 
outer core of a neutron star reduces to the search for 
the function Ejs[{nn,np). A number of ways to address 
this problem have been proposed based on a variety of 
theoretical physics methods, viz, the Brueckner-Bethe- 
Goldstone theory, the Green's function method, varia- 
tional methods, the relativistic mean field theory, and 
the density functional method 25, §5.9]. The Akmal- 
Pandharipande-Ravenhall (APR) model known in sev- 
eral varia nts is currently regarded as the most reliable 
one [ll5j . The APR model uses the variational princi- 
ple of quantum mechanics, under which an energy min- 
imum for the trial wave function is sought. This func- 
tion is constructed by applying the linear combination 
of operators describing admissible symmetry transfor- 
mations in the coordinate, spin, and isospin spaces to 
the Slater determinant consisting of wave functions for 
free nucleons. APR variants differ in the effective po- 
tentials of nucleon-nucleon interaction used to calculate 
the mean energy. The potentials borrowed by the au- 
thors from earlier publications take the modern nuclear 
theory into account and their parameters are optimized 
so as to most accurately reproduce the results of nuclear 
physics experiments. We note that the addition of an 
effective three-particle nucleon-nucleon potential to the 
two-particle one ensures a remarkably close agreement 
between theory and experiment. 

The effective functional of nuclear matter energy den- 
sity was used to construct another known equation of 
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Figure 2: FPS, Sly, and APR equations of state for the 
core of a neutron star. Bold dots on the curves corre- 
spond to the maximahy possible density in a stationary 
star. 



state, SLy |116j . Calculations based on the SLy equa- 
tion are less detailed but easier to use than in the APR 
model. This equation is constructed in accordance with 
the s ame scheme as the well-known FPS equation of state 
[117| . which was especially popular in the 1990s in calcu- 
lations of the astrophysical properties of neutron stars. 
The main difference between SLy and FPS lies in the 
specification of parameters of the effective energy density 
functional accounting for current experimental data. An 
important advantage of both models over many others is 
their applicability not only to the stellar core but also to 
the crust, which allows determining the positio n of the 
crust-core interface in a self-consistent manner |118j . 

Note by the way that there are a conve nient 
parametrization for the APR equation of st ate |119l | 
and explicit fitting expressions for Sly jl20| for the 
dependence of pressure on density and the so-called 
pseudoenthalpy, a convenient parameter for cal culat ing 
the properties of rapidly rotating neutron stars [121j . 

Figure shows P{p) dependences for models FPS, 
SLy, and APR. Comparison of FPS and Sly shows that 
a more exact account of current experimental data makes 
the P{p) dependence steeper and the equation of state 
stiffer. Bold dots correspond to the density in the center 
of a neutron star with M = Mmax for each of these 
equations; the segments of the curves to the right of 
these dots cannot be realized in a static star. 

A common drawback of the above models is the appli- 
cation of a Lorcntz non-invariant theory to the descrip- 
tion of hadrons. Such a description becomes a priori in- 



correct in the central part of the core, where the speeds of 
nucleons on the Fermi surface may constitute an appre- 
ciable fraction of the speed of light. The same drawback 
is inherent in all other aforementioned approaches, with 
the exception of the relativistic mean field theory. This 
theory, suggested in the 1950s, was especially popular in 
the 1970s [122]. It has a number of appealing features. 
Specifically, its Lorentz invariance guarantees the fulfill- 
ment of the condition that the speed of sound does not 
exceed the speed of light, which is subject to violation in 
some other models. But the assumption of spatial uni- 
formity of meson field sources underlying this theory is 
valid only if rih ^ 100 no [1^. The matter density nec- 
essary for this condition to be satisfied is much higher 
than that in the interior of a neutron star. Therefore, 
the equations of state in the core of neutron stars are cal- 
culated realistically based on nuclear interaction models 
that are not Lorentz invariant but are still applicable to 
the largest part of the stellar core. 

Along with the P{p) dependence, it is important to 
know relative abundances of various particles. In par- 
ticular, the dependence of the proton fraction Xp in the 
neutron-proton-electron-muon (npefi) matter on density 
p. The fact is that the principal mechanism behind neu- 
trino energy losses in the outer core of a neutron star is 
the so-called modified Urea processes (in short, Murca, 
after K P Levenfish) consisting of consecutive reactions 
n + N p+N + e + Pe and p+ N + e ^ n + N + i/e, where 
iV = n or A'^^p is a nucleon- mediator ("an active spec- 
tator", according to Chiu & Salpeter (l^)- The involve- 
ment of the mediator distinguishes the Murca processes 
from ordinary processes of beta-decay and beta-capture, 
referred to as direct Urea processeso- If Xp < Xc, where 
the value of Xc varies from 0.111 to 0.148 depending 
on the muon abundance, then the energy and momen- 
tum conservation laws cannot be fulfilled simultaneously 
without the participation of a mediator nucleon in the 
Urea process, keeping in mind that the momenta of the 
involved strongly degenerate neutrons n, pr otons p, and 
electrons e lie near their Fermi surfaces |125j . If Xp > Xc, 
then direct Urea processes much more powerful than 
Murca come into play. For this reason, the neutron star 
suffers enhanced cooling if Xp exceeds Xc- 

In different variants of the APR model, the fraction of 
Xp grows nonmonotonically from « 0.01 - 0.02 at Ub — uq 
to Xp « 0.16-0.18 at Ub > 1.2 fm^^ Therefore, ii the 
APR model holds, the density of matter in the center 
of a massive (M > 1.8 Mq) neutron star is such that 
it switches on direct Urea processes and speed up its 
cooling. In contrast, direct Urea processes are impossible 



^T he term Urea pro cess was coined by Gamow and Schonberg 
[123l |. Gamow recalled [124| : "We called it the Urea process par- 
tially to commemorate the casino in which we first met and par- 
tially because the Urea process results in a rapid disappearance 
of thermal energy from the interior of a star similar to the rapid 
disappearance of money from the pockets of the gamblers in the 
Casino da Urea. Sending our article on the Urea process for publi- 
cation in the Physical Review, I was worried that the editors would 
ask why we called the process 'Urea'. After much thought I de- 
cided to say that this is the short for unrecordable cooling agent, 
but they never asked." 
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Figure 3: Equations of state for the core of a neutron star 
according to the variants of the model proposed in [127l | 
with (BGNIHI and BGN1H2) and without (BGNI) hy- 
perons. Bold dots on the curves indicate the maximally 
possible density in a stationary star. The data and no- 
tations are borrowed from [2511 . 



(for stable stars) in the Sly model. 

4.2 The inner core and hyperons 

Strong gravitational compression in the interior of a neu- 
tron star is likely to provoke the conversion of nuclcons 
into hyperons if such conversion may reduce the energy 
density at a given ni,. The process is mediated by the 
weak interaction with a change of strangeness (quark 
flavor). According to modern theoretical models, the 
conversion is possible at p > 2po. 

The equation of state containing hyperons is calcu- 
lated as described for the npefj, case in Sect. 14. 1[ but the 
equations for the chemical potentials are supplemented 
with new ones for the equilibrium conditions with respect 
to the weak interactions. The lightest baryons make an 
octet of two nucleons {p and n with zero strangeness 
S = 0), four hyperons with S = -1 (A°, S", E", and 
S+), and two hyperons with S — ~2 (S° and S~). 
Here, they are listed in the order of increasing mass. 
Under normal conditions, hyperons decay for fractions 
of nanoseconds. But in matter composed of degenerate 
neutrons, /i„ increases with increasing density. When /i„ 
reaches a minimal chemical potential of a hyperon given 
by its mass, this hyperon becomes stable because the de- 
cay reaction ceases to be thermodynamically favorable. 

Some clarification is needed here. Although A° is the 
lightest of all hyperons, S" is the first to be stabilized 



as the density increases. It occurs so, because the decay 
of the hyperon yields an electron, in conformity with 
the equilibrium condition /zs- = /x„-|-yUe. The electrons, 
like neutrons, are strongly degenerate, and their chemi- 
cal potential /ig is equal to the Fermi energy; the addition 
of this energy to fin permits the equilibrium condition to 
be satisfied at a smaller density (as was first noticed by 
Salpeter in 1960 (126] ). Similarly, the necessity of sub- 
tracting /Ze from can make formation of S"*" hyperons 
disadvantageous. However, electrons are gradually re- 
placed by S~ hyperons, and fie decreases as the density 
increases. Due to this, many theoretical models predict 
the appearance of S"*" hyperons at a sufficiently high den- 
sity, rife > 5 no- In the general case, both electrons and 
muons are gradually substituted by negatively charged 
hyperons as the density increases. In the models predict- 
ing a high hyperon concentration, leptons disappear at 
ri}, > 1 im~^ and the so-called "baryonic soup" is cooked 
with high average strangeness per baryon (almost —1 in 
the central parts of maximum-mass stars). 

The current theory lacks a rigorous description of 
nuclcon-hyperon and hyperon-hyperon interactions. 
This uncertainty is aggravated by the uncertainty 
arising from the choice of the mode of description of 
multiparticle interactions; this results in a great variety 
of model equations of state for the inner core of a 
neutron star. Figure [3] exe mplifies three of the equations 
of state proposed in [127f . The solid curve corresponds 
to the so-called minimal model disregarding hyperons. 
The dashed and dashed-dotted curves correspond to 
two models with hyperons. If the density exceeds 
the threshold for the appearance of new particles, the 
equation of state is noticeably softened, as is natural 
when part of the strongly degenerate high-energy neu- 
trons are replaced by slow heavy hyperons. However, 
the magnitude of the effect depends on the details of 
hyperon-nucleon and hyperon-hyperon interactions. 

4.3 Phase transformations and decon- 
finement 

As the density p increases above the nuclear density pa 
matter may undergo phase transitions to qualitatively 
new states regarded as exotic from the standpoint of ter- 
restrial nuclear physics; the very existence of these states 
depends on the concrete features of strong interactions 
and the quark structure of baryons. 

4.3.1 Meson condensation 

It has been known since the mid-1960s [tH that the core 
of a neutron star must contain it mesons (pions), i.e., the 
lightest mesons. Bose condensation of pions in nuclear 
matter is usually hampered by strong pion-nucleon re- 
pulsion. However, it was shown in [128 131J that collec- 
tive excitations (pion-like quasiparticles) may arise in a 
superdense medium and condense with the loss of trans- 
lational invariance. Further studies revealed the possi- 
bility of creating different phases of the pion condensate 
and the importance of correlations between nucleons for 
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its existence. It was shown that short-range correlations 
and the formation of ordered structur es in the dense mat- 
ter interfere with pion condensation [132j |. 

Kaons (i4r-mesons)are the hghtest strange mesons. 
They appear in the core of a neutron star as a resuh 
of the processes e + N ^ + N + and n + — > 
p -t- K" -t- N , where iV is a nucleon whose participation 
ensures the momentum and energy conservation in the 
degenerate matter. The possibihty of Bose condensa- 
tion of kaons at p > 3po first understood in the 1980s 
133 has thereafter been studied by many authors (see 
134| for a review). In a neutron star, it involves K~- 
like particles, by analogy with pion condensation. These 
particles have a smaller mass than isolated K mesons; 
it is this property that makes their Bose condensation 
possible. A method for the theoretical description of 
a kaon condensate taking the effects of strong interac- 
tion in baryonic matter into account was developed in 
[135j . Formation of the kaon condensate depends on the 
presence of hyperons and strongly affects the properties 
of the nucleon component of matter. Kaon condensa- 
tion, like pion condensation, is accompanied by the loss 
of translational invariance. The condensate forms via 
first- and second-order phase transitions, depending on 
the stren gth o f the force of attraction between kaons and 
nucleons [136j | . Both pion and kaon condensations make 
the equation of state much softer. 

4.3.2 Quark deconfinement 

Because hadrons are made of quarks, the fundamental 
description of dense matter must take the quark degrees 
of freedom into account. Quarks cannot be observed 
in a free state when their density is low because they 
are held together (confine ment ) by the binding forces 
enhanced at low energies |137l |. As the density (and 
hence, the characteristic energies of the particles) grows, 
baryons fuse to form quark matter. In 1965, Ivanenko 
and Kurdgelaidzc 138] suggested that neutron stars have 
quark cores. With the advent of quantum chromody- 
namics, calculations of quark matter properties were per- 
formed in terms of the perturbation theory using the 
nonintera cting quark model as the initial approximation 
[139l ll4Cll |. However, the use of this theory is limited to 
energies ^ 1 GeV, while the chemical potential of par- 
ticles in neutron stars does not reach such high values. 
More quark matter models were proposed, and the su- 
perfluidity of this matter associated with quark Cooper 
pairing was considered (see the references in pSl §7.5]). 
These models were used to explore quark stars and hy- 
brid stars, i.e., neutron stars with cores mad e of quark 
matter 0. For instance, the authors of |141| predicted 
a series of phase transitions in the interior of a hybrid 
star with sequential deconfinement of quark flavors at 
Ub ~ 0.25, 0.5-0.8, and 1.1-1.8 fm'^. 

All published models of phase transitions in the cores 
of neutron stars have serious drawbacks. The quark and 
baryonic phases are typically treated in the framework of 
different models and cannot therefore be described self- 
consistently. Calculations from the perturbation theory 



are unrealistic at those relatively small densities at which 
they predict a phase transition. For this reason, the ex- 
istence of a quark core in neutron stars cannot be proved 
theoretically. However, it may be hoped that this goal 
will be achieved based on the analysis of observations of 
compact stars. 

4.3.3 Mixed phases 

First-order phase transitions can be realized via a state 
in which one phase co-exists with another in the form 
of dr oplets. Such phase transitions, called noncongruent 
|142| . have been cons ider ed in connection with compact 
stars since the 1990s jl43j .The coexistence of two phases 
in the core of a neutron star is possible thanks to the 
abandonment of the implicit assumption of electroneu- 
trality of each individual phase. In the mixed phase, 
the electric charge of one component is on the average 
compensated by that of the other, and the matter struc- 
ture is determined by the balance of surface tension at 
the boundaries between droplets, the energy density of 
baryonic matter, the kinetic energy of the constituent 
particles, and electrostatic energy. Mixed states are fea- 
sible for both meson condensation and baryon dissocia- 
tion into quarks. 

4.3.4 Crystalline core 

The early models of neutron stars assumed that strong 
short-range neutron-neutron repulsion results in th e for - 
mation of the solid inner core of a neutron star [144l |. 
as mentioned in review by Ginzburg fT9l |. In subse- 
quent works, it was taken into account that the nucleon- 
nucleon interaction occurs by an exchange of vector 
mesons, which gives rise to the effective Yukawa poten- 
tial. As the calculations became more exact by the late 
1970s, it was understood that the realistic effective po- 
tentials of neu tron- neutron interactions do not lead to 
crystallization |145j . 

An alternative possibility of crystallization arises from 
the tensor comp onent of the mid-r ange nucleon- nucleon 
attraction [146l |. It was shown in |147| that tensor in- 
teraction can lead to structures in which neutrons are 
located in a plane with oppositely oriented spins, each 
such plane hosting oppositely directed proton and neu- 
tron spins - so-called alternating spin (ALS) structures. 
If the gain in the binding energy during formation of 
an ALS structure exceeds the loss of the kinetic energy 
of the particles, then the structure may become ener- 
getically favorable and a phase transition into this state 
occurs. 

Moreover, given a low enough abundance of protons in 
baryonic matter, Xp < 0.05, their localization may occur , 
accompanied by modulation of the neutron density jl48{ . 
Under certain conditions, mixed phases may also be just 
ordered into periodic structures (see monograph and 
the references therein). 

In other words, there are numerous hypotheses re- 
garding the structure and composition of the neutron 
star core, differing in details of microscopic interactions 



14 



A Y POTEKHIN 



and theoretical models for their description. We consider 
how they treat the parameters of neutron stars. 

4.4 Relation to observations 

Manifestations of the properties of a neutron star core 
can be arbitrarily categorized into dynamic and qua- 
sistatic. The former are due to relatively fast processes 
inside the star. For example, a phase transition in the in- 
ner core may occur not only at the star birth but also as 
it goes through its evolution, e.g., during cooling (when 
the temperature decreases to below the critical one) or 
rotation slowdown (when the central pressure increases 
due to a reduction in centrifugal forces). Such a phase 
transition leads to a starquake with the release of ther- 
mal energy, a burst of neutrinos, excitation of crustal 
oscillations, and an abrupt change in the rotational ve- 
locity due to the altered moment of inertia [54| . All these 
effects could be possible to record and mea sure under 
favorable conditions. The authors of jl49| attributed 
sharp jumps of pulsar rotational periods (glitches) to 
such phase transitions. Starquakes and glitches may be 
a consequence of occasional adjustment of the crust ro- 
tational velocity to the rotation rate of the superfluid 
component of the nucleon liquid [H, HI, [111 . 

Quasistatic manifestations include, inter alia, the ef- 
fects of the core structure and the physical properties of 
its superdense matter on the theoretical radius and cool- 
ing rate of the star. The influence on the cooling, and 
hence the effective surface temperature, originates from 
the difference in the rates and mechanisms of neutrino 
losses in individual core models ^26*, '3^ . The influence 
on the relation between the stellar radius R and mass 
M is realized via the function P{p), in accordance with 
the TOV equation. Figure |4] exemplifies the R{M) de- 
pendence for six equations of state of the neutron stars 
shown in Figs. [2] and |3l and one equation for a quark 
(strange) star (with the use of data from 25] ) . It can be 
seen that quark stars must have smaller masses and radii 
than typical neutron stars. Solid dots at the ends of the 
curves correspond to the maximum mass of a stationary 
star for each equation of state. If a neutron star of a 
higher mass is discovered, the equation in question can 
be discarded. 

Besides the equations of state, there are general theo- 
retical constraints on the possible values of masses and 
radii. Evidently, the radius R of any star must not be 
smaller than r^; otherwise, we are dealing with a black 
hole. Moreover, it can be shown [25[ that the condition 
Vs < c, where Vs is the speed of sound in a local refer- 
ence frame and c is the speed of light in the vacuum, im- 
posed by the special theory of relativity and the causal- 
ity principle, requires that the relation R > 1.412 r^, be 
satisfied, which excludes the point (M, R) from entering 
inside the hatched triangle in Fig. |4l 

An additional limitation is needed for the gravitation 
of a rotating star to overcome the centrifugal accelera- 
tion. Clearly, the radius must not be too large. In Fig. HI 
the largest values of the radius in the dependence on the 
mass M at a given rotation period P are shown by al- 



ls 



B 

A 10 
5 



■ 

r 
1 


1 1 1 1 1 1 1 1 1 1 1 


1 ,1 

CO 

+ 
n 
o 
ai 


/ ^- 




1 

- 




^ a: 
m 
P-, 

/ 


'r 


■ 






'V\ — 






1: 
2: 
3: 


APR / 
SLy 

FPS / 


~ I 
1 . 


V 3 2 • * 1^ 






4: BGNl ,,<<>^> 
5: BGNIHI ,<<$$^^^ 
6: BGN1H|<<^^J?^;^^ 








- / 











0.5 1 1.5 2 2.5 



M/M, 



Figure 4: Mass dependence of the compact star radius. 
Curves 1-6 correspond to the equations of state of neu- 
tron stars shown in Figs. [2] and [H The dotted curve 
corresponds to one of the feasible equations of state of 
a quark star. The hatched triangular area is forbidden 
by the causality principle. The crosshatched triangu- 
lar area lies below the event horizon. Three short- and 
long-dashed curves show additional constraints for ro- 
tating neutron stars: the area under the corresponding 
curve is permitted at the rotation periods indicated (1.4, 
1, and 0.5 ms). The straight line R = Sr^ corresponds 
to the minimal stability radius of the circular orbit of 
a test particle around the neutron star of a given mass. 
The wide vertical strip encompasses the masses of binary 
neutron stars measured with an error less than 0.1 M© at 
confidence level 2a [25| : the narrow vertical strip corre- 
sponds to the mass of the PSR J1903-I-0327 millisecond 
pulsar 02] ■ (See the Note added in proof.) 



ternating short and long dashes for P = I ms, 1.4 ms, 
and 0.5 ms. We see that the period P = 1.4 ms (the 
shortest of the periods observed to date) does not place 
any serious constrains on R. On the other hand, the pul- 
sar period 0.5 ms is incompatible with any of the known 
theoretical equations of state of dense matter (the de- 
tection of such a period in the radiation fro m th e 1987A 
supernova remnant was reported ii i 19 89 |150| . but it 
later proved to be a technical error [151| ). 

The wide vertical strip in Fig. |4] depicts the range of 
exactly measured masses of neutron stars in binary sys- 
tems made up of a pulsar and another neutron star The 
narrow vertical strip corresponds to the estimated mass 
of PSR J1903-H327 mentioned in Sect. O Confirmation 
of this estimate would make the choice between theoreti- 
cal models much more definitive. For example, it follows 
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from the figure that the existence of a star of such mass 
implies the absence of hyperons in the model [127]. 

Were R and M known exactly for a certain compact 
star, it would probably permit choosing one of the equa- 
tions of state as the most realistic one. Unfortunately, 
the current accuracy of measurement of neutron star 
radii leaves much to be desired. Determination of stel- 
lar masses and radii requires a reliable theoretical de- 
scription of the envelopes that influence the star surface 
temperature and the formation of the emitted radiation 
spectrum. We return to this issue in Sect. 17.51 

5 Envelopes 

Envelopes of a neutron star are divisible into the solid 
crust in which atomic nuclei are arranged into a crystal 
and the liquid ocean composed of the Coulomb fluid. 
The crust is subdivided into the inner and outer parts. 
In the former, the nuclei are embedded in a sea of free 
neutrons and electrons, while the latter contains no free 
neutrons. A neutron star can have a gaseous plasma 
atmosphere at the surface, while the stellar core may 
be surrounded by a liquid crystal mantle topped by the 
crust. 

5.1 Inner crust 

The inner crust s normally ~ 1 - 2 km thick. Its density 
increases from pdrip ~ (4-6) x 10^^ g cm~^, , at which 
neutrons begin to "drip" from the nuclei, to ^ 0.5pO) 
when the atomic nuclei fuse into a homogeneous mass. 
The nuclear chemical equilibrium with respect to beta- 
capture and beta-decay reactions in the inner crust ac- 
counts for the matter composition that cannot be repro- 
duced under laboratory conditions (neutron-rich heavy 
nuclei embedded in a fluid composed of neutrons and 
electrons). The physics of such matter is fairly well de- 
scribed in [152 ]. The neutrons in a large portion of the 
inner crust are superfluid; according to theoretical esti- 
mates, the critical superfluidity temperature varies with 
density and reaches billions of degrees or an order of 
magnitude higher than the typical kinetic temperature 
of matter in the inner crust of a neutron star. 

The pressure in the inner crust of a neutron star is 
largely created by degenerate neutrons. However, super- 
fluidity may decrease their heat capacity and is therefore 
responsible for the decisive contribution of atomic nuclei 
to the thermal capacity of the inner crust. The nuclei 
make up a crystal lattice, formed essentially by Coulomb 
interaction forces (Coulomb or Wigner crystal). An ade- 
quate description of their contribution is possible by con- 
sidering collective vibrational excitations (phonon gas). 
Because the electrons are relativistic and highly degen- 
erate particles, their contribution to the heat capacity 
of the inner crust is insignificant unless the temperature 
is too low. However, it may appreciably increase when 
the temperature of the Coulomb crystal decreases much 
below the Debye temperature, imposing a "freeze-out" 
on phonon excitations §2.4.6]. 



The electric conductivity in the inner crust is due 
largely to electrons, whose scattering plays an impor- 
tant role. The scattering on phonons of a crystal lattice 
prevails at relatively high temperatures, and that on lat- 
tice defects or admixtures is responsible for the residual 
resistance at low temperatures. Ions (atomic nuclei) in- 
corporated into the crystal lattice make no appreciable 
contribution to conductance. At the same time, thermal 
conductivity is due to phonons and neutrons, besides 
electrons whose scattering is governed by the same mech- 
anisms that operate in the case of electric conductivity 
supplemented by electron-electron collisions. Phonons 
may become the main heat transfer agents in the pres- 
ence of lattice defects and admixtures hampering the 
participation of electrons in this process [1531] . Neutrons, 
especially superflu id on es, may also serve as heat carriers 
in the inner crust [154| . 

5.2 Mantle 

The core of a neutron star may be separated from the 
bottom of its inner crust by a layer that contains ex- 
otic atomic nuclei and is called the mantle 155]. In the 
liquid-drop model, the spherical shape of the atomic nu- 
cleus is energetically advantageous at a low density, since 
it minimizes the surface energy. However, the contribu- 
tion from the Coulomb energy at a higher density may 
change the situation. The mantle consists of a few layers 
containing such phases of matter in which atomic nuclei 
are shaped not like spheres but rather likes cylinders 
(the so-called "spaghetti" phase), plane-parallel plates 
( "lasagna" phase) , or "inverse" phases composed of nu- 
clear matter with entrapped neutron cylind ers ("tubu- 
lar" phase) and balls ( "Swiss cheese" phase) fl52] . Such 
structures for c ollap sing cores of supernovae were firs t 
conjectured in |156j . and for neutron stars in |157f . 
Whereas the spherical nuclei constitute a 3D crystalline 
lattice, the mantle has the properties of a liquid crystal 
[155]. Direct Urea processes of neutrino emission can be 
allowed in the mantle [158], while they are unfeasible in 
other stellar envelopes; their high intensity can enhance 
the cooling of the neutron star. 

Not all current equations of state of nuclear matter 
predict the mantle; some of them treat such a state as 
energetically unfavorable. The mantle hypothesis ap- 
pears in the FPS model, but not in the most modern Sly 
model. 

5.3 Outer crust and its melting 

The outer shells of a neutron star are hundreds of me- 
ters thick and consist of an electron-ion plasma that is 
completely ionized, that is, consists of ions in the form 
of atomic nuclei and strongly degenerate free electrons 
(probably except a several-meter-thick outer layer with 
the density below 10^ g cm~^). Then the total pressure is 
determined by the pressure of degenerate electrons. The 
electrons become relativistic (with the Fermi momentum 
Pf comparable to WeC, where rrie is the electron mass) 
at p > 10^ g cm~^ and ultrarelativistic {pp ^ mec) at 
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Figure 5: Schematic cross section of neutron star envelopes. From bottom up: core, mantle, inner crust, outer 
crust, ocean, atmosphere. Right-hand side: composition of these layers; left-hand side: characteristic values of 
density logarithm and depth from the surface. 



p ^ 10^ g cm~^. At such densities, ions give rise either 
to a Coulomb liquid (whose properties mostly depend 
on Coulomb interactions between ions) or to a Coulomb 
crystal. 

The electron Fermi energy in deep-lying layers of the 
outer envelopes increases so as to enrich nuclei with neu- 
trons by virtue of beta-captures. Finally, the inner-outer 
crust interface forms at p — pdrip, where free neutrons 
appear. 

The external boundary of the outer crust normally 
coincides with the crystallization point of the Coulomb 
liquid making up the neutron star ocean. The posi- 
tion of this point is given by the density dependence 
of the Coulomb crystal melting temperature. In the so- 
called one-component Coulomb plasma model disregard- 
ing electron-ion interactions and treating ions as classical 
point particles, formation of the Coulomb crystal is de- 
fined by the equality F = 175 or (in a more realistic ap- 
proach) r ~ 100-200 dH. Here, F = (Z e)^ / {akBT)is 
the Coulomb coupling parameter characterizing the re- 
lation between the potential Coulomb energy of the ions 
and their kinetic energy, a — (47r?ij/3)^^/'^ is the ion- 
sphere radius, Ui is the ion volume density, and is the 
Boltzmann constant. The melting point for a typical 
neutron star envelope hes at pm ~ 10^ — 10^ g cm~^, de- 
pending on its thermal structure (i.e., temperature vari- 
ations with depth related to the star age and past his- 
tory). However, a cold enough neutron star may lose 
both the atmosphere and the ocean in a superstrong 
magnetic field; in such a case, the external boun dary 
of the crust coincides with the stellar surface (see |160l | 
for discussion and references). 

5.4 Ocean 

The bottom of the neutron star ocean is located at the 
melting point with the density pm, while its surface 
is arbitrary because of the lack of a clear-cut ocean- 
atmosphere interface on a typical star. An exception, 
as in the case of the solid crust, is neutron stars having 



a rather strong magnetic field, which may be responsible 
for the absence of an optically thick atmosphere and its 
substitution by a liquid boundary. Most of the ocean 
consists of atomic nuclei surrounded by degenerate elec- 
trons. Therefore, in the general case, we speak of ions 
surrounded by electrons, with the understanding that 
ions mean both completely and partially ionized atoms. 

The ocean matter is the Coulomb liquid, most of which 
is strongly coupled, i.e., F 3> 1. One of the main prob- 
lems in theoretical studies of such matter is adequate 
consideration of the influence of microscopic correlations 
between ion positions on the macroscopic physical char- 
acteristics of the matter being investigated, such as equa- 
tions of state 159] and kinetic coefficients [161j . 

5.5 Atmosphere 

The stellar atmosphere is a layer of plasma in which the 
thermal electromagnetic radiation spectrum is formed. 
The spectrum contains valuable information about the 
effective surface temperature, gravitational acceleration, 
chemical composition, magnetic field strength and ge- 
ometry, and mass and radius of the star. The geometric 
thickness of the atmosphere varies from a few millimeters 
in relatively cold neutron stars (effective surface temper- 
ature Toff ^ 10^-^ K) to tens of centimeters in rather hot 
ones (Toff ^ 10^'^ K). In most cases, the density of the 
atmosphere gradually (without a jump) increases with 
depth; however, as mentioned above, stars with a very 
low effective temperature or superstrong magnetic field 
have either a solid or a liquid condensed surface. 

The deepest layers of the atmosphere (its "bottom" 
being defined as a layer with the optical thickness close 
to unity for the majority of outgoing rays) may have 
the density p from ^ 10~^ to ~ 10^ g cm~^, depending 
on the magnetic field B, temperature T, gravity g, and 
the chemical composition of the surface. The presence 
in the atmosphere of atoms, molecules, and ions having 
bound states substantially alters absorption coefficients 
of electromagnetic radiation and thereby the observed 
spectrum. 
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Although the neutron star atmosphere has been in- 
vestigated by many researchers for several decades, these 
studies (especially concerning strong magnetic fields and 
incomplete ionization) are far from being completed. 
For magnetic fields B ~ 10^^ - 10^'' G, this problem is 
practically sol ved only for hydrogen atmospheres with 
T^cff > lO^-S K [162', Tes]. The bound on Tcff from below 
is related to the requirement of smallness of the contri- 
bution from molecules compared to that from atoms, the 
quantum mechanical properties of molecules in a strong 
magnetic field being poorly known. For B ^ 10^^- 
10" G and 10^-^ K < Tcff < 10^ K, there are models 
of partly ionized atmospheres composed of carbon, oxy- 
gen, and nitrogen [164]. Here, the restriction on Tefi 
from above arises from the rough interpretation of ion 
motion effects across the magnetic field that holds at 
low thermal velocities (see Sect. 16. 3p . 

6 Magnetic fields 

6.1 Magnetic field strength and evolu- 
tion 

As discussed in the Introduction, the majority of 
currently known neutron stars have magnetic fields 
unattainable in terrestrial laboratories, with typical 
values B ~ 10® - 10^^ G at the surface, depending on 
the star type. The field strength inside a star can be 
even higher. For example, certain researchers propose 
explaining the energetics of AXP and SGR in terms of 
core magnetic fields as high as B ~ 10^^ - 10^^ G at the 
birth of the neutron star (see |165l | and the references 
therein). The theoretical upper bound obtained numer- 
ically in |121| is consist ent with the estimate from the 
virial theorem [HI, [Till: max(B) - 10^® G. 

A number of theoretical models of field generation 
have been proposed suggesting the participation of dif- 
ferential rotation, convection, magneto-rotational insta- 
bility, and thermomagnetic effects either associated with 
supernova explosion and collapse or occurring in young 
neutron stars (s ee Il68| '). . Specifically, the "a-0- 
dynamo model" [161 Il70l| assumes that the core of a 
neutron star born with a sufficiently short (millisecond) 
rotation period acquires a toroidal magnetic field up to 
B ~ 10^^ G due to differential rotation, while the pulsar 
magnetic field is generated by means of convection at 
the initial rotation periods > 30 ms. However, none of 
the proposed models is able to account for the totality 
of currently available neutron star data. 

Electric currents maintaining the stellar magnetic field 
with the involvement of differential rotation circulate ei- 
ther in the inner crust or in the core of a neutron star, 
i.e., where electric conductivity is high enough to pre- 
vent field decay for a time comparable with the^age of 
known pulsars. It was shown as early as 1969 (l7H that 
the characteristic time of Ohmic decay of the core mag- 
netic field may exceed the age of the Universe. For a 
magnetic field originating in the core of a neutron star, 
proton superconductivity stipulates its existence in the 



form of quantized magnetic tubes (Abrikosov vortices, 
or fiuxoids) having a microscopic transverse dimension. 

The magnetic field of a neutron star changes in the 
course of its evolution, depending on man y fac tors and 
inter- related physical processes (see, e.g., |172j and the 
references therein). Specifically, the field undergoes 
ohmic decay and a change in configuration under the 
effect of the Hall drift; also, magnetic force lines recon- 
nect during starquakes. Thermoelectric effects, as well as 
the dependence of components of the thermal and elec- 
trical conductivity tensors, and plasma thermoelectric 
coefficients on temperature and the magnetic field, are 
responsible for the interrelation between magnetic and 
thermal evolution [H, . Accretion ca n als o strongly 
affect the near-surface magnetic field [2l|, |172| . 

The evolution of a magnetic field generated by 
Abrikosov vortices is to a large extent dependent on 
their interaction with other core components, such as 
Feynman-Onsager vortices in the neutron superfluid 

[1, Il71| . and conditions at the core boundary^ i.e., 

interactions of these vortices with crustal matter [S^l . 

6.2 Landau quantization 

The motion of free electrons p erpe ndicular to the field is 
quantized into Landau levels |173j . Their characteristic 
transverse scale is the magnetic length a™ = (hc/eBy/^ , 
and the inter-level distance in a nonrelativistic theory 
is the cyclotron energy huJc = 11.577 B12 keV, where 
cOc = eB / (m,f,c) is the electron cyclotron frequency (the 
notation B12 — i?/(10^^ G)is introduced here). The di- 
mensionless parameters characterizing a magnetic field 
in relativistic units b and atomic units 7 are 

b = ;ic^,/(mec2) = Bi2/44.14, (7) 

where ob is the Bohr radius. We call a magnetic field 
strong if 7 ^ 1 and superstrong if 6 > 1. In the 
relativistic theory, the energies of Landau levels are 
En = rriec^ (Vl + 267V - 1) {N = 0, 1, 2, . . .). 

In a superstrong field, specific effects of quantum elec- 
trodynamics, such as electron-positron vacuum polariza- 
tion in an electromagnetic wave field, become significant. 
As a result, the vacuum acquires the properties of a bire- 
fringent medium, which, at 6 > 1, markedly affects the 
radiation spectrum formed in the atmosphere of a neu- 
tron star 174, 175]. 

For ions with charge Ze and mass rrii = Aruu, where 
m„ = 1.66 X 10~^* g is the unified atomic mass unit, 
the cyclotron frequency equals ujd ~ \Ze\B / [mic), the 
cyclotron energy hwd = 6.35 [Z/A) B12 eV, and the pa- 
rameter that characterizes the role of relativity effects 
bi = hujci/imic^) = 0.68 x 10"® {Z/A'^)Bi2. The small- 
ness of bi allows one to ignore the relativistic effects for 
ions in the atmosphere of a neutron star. 

The motion of electrons along a circular orbit in a 
magnetic field in the classical theory leads to cyclotron 
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radiation at the frequency Wc and the formation of a cy- 
clotron hne in the magnetic atmosphere. In quantum 
theory, the cyclotron line corresponds to transitions be- 
tween the adjacent Landau levels. Transitions between 
distant Landau levels result in the formation of cyclotron 
harmonics with energies E]\[. The discovery of the first 
cyclotron line with the energy 58 keV in the spectru m of 
the X-ray pulsar in the binary Hercules X-1 in 1978 |176j 
gave a stunning argument in favor of the idea of neutron 
star magnetic fields. A number of similar systems with 
synchrotron lines are known presently. The spectra of 
several X-ray pulsars in binary systems exhibit cyclotron 
harmonics with N ^ 1 (see 177, 178]); th e obs ervation 
of up to four harmonics has been reported [179j |. 

The spectra of isolated neutron stars may likewise 
contain cyclotron lines. Electron cyclotron lines can be 
seen in the thermal spectral range from 0.1 to 1 keV at 
B - 10^°- 10" G and ion cyclotron lines at B ~ lO^^- 
10^"* G. There is a hypothesis that absorption lines in the 
spectrum of CCO IE 1207.4-5209 can be explained by 
the cyclotron mechanism |180l Il81j . It should be noted 
that ion cyclotron harmonics, unlike electron ones, are 
too weak to be observed [9^. 

The effect of Landau quantization on plasma proper- 
ties is significant when the cyclotron energy is not too 
small compared with the thermal (ksT) and Fermi (ep) 
energies. If huc is much higher than these two ener- 
gies, most electrons in thermodynamic equilibrium are 
at the ground Landau level. In this case, the field 
is called strongly quantizing. If, at contrast, k^T or 
epis much greater than the energy difference between 
the adjacent Landau levels, the field is nonquantiz- 
ing. The smallness condition on the thermal energy 
compared with the cyclotron one can be written as 
tujJc/{k^T) « IMB12/TQ > 1 for the electrons and 
huci/kBT « 0.0737 {A/Z)Br2/Te > 1 for the ions. The 
second condition (higher huc than ep) imposes a bound 
on density. Degenerate electrons occur at the ground 
Landau level if their number density Ue is smaller than 
ub = (tt^v^ a^)~^. Therefore, the field will be strongly 
quantizing at p < ps, where 

PB = miUB/Z ^7x10^ (A/Z) Sff g cm'^. (9) 

At p > pb the field is weakly quantizing, and at p 3> ps 
it can be treated as nonquantizing. Estimates, analogous 
to (HI), are feasible also for other fermions [23 §5.17]. 

We note that a nonquantizing magnetic field has no 
effect on the equation of state (Bohr - van Leeuwen the- 
orem) . It follows from Fig. [S] with account of ^ that 
magnetic fields inherent in neutron stars are strongly 
quantizing in the atmosphere and can be quantizing in 
both the ocean and the outer crust; however, even at 
limit values B ~ 10^^ G, the magnetic fields do not af- 
fect the stellar core equation of state. These conclusions 
are based on simple estimates, but they are confirmed 
by calculations of the nuclear matter equation of state 
in superstrong magnetic fields |182j . 



6.3 Atoms and ions in magnetic atmo- 
spheres 

The atmosphere of a neutron star contains atoms, 
molecules, and atomic and molecular ions having bound 
states. Strong magnetic fields markedly a ffect their 
qua ntum mechanical properties (see reviews [Tea, EH, 
|184| ) It was suggested soon after the discovery of pul- 
sars |185| . that at equal temperatures, there should be 
more atoms in the neutron star atmosphere at 7 ^ 1 
than at 7 < 1, because in a strong magnetic field, 
the binding energies of their ground state and a cer- 
tain class of excited states (so-called tightly bound states) 
markedly increase and the quantum mechanical size de- 
creases. For instance, the ground-state energy of an 
H atom at _B ^ 10^^-10^'' G can be roughly esti- 
mated as E ^ 200(lni?i2)^ eV. In all the states at 
7^1, the electron cloud acquires the form of an ex- 
tended ellipsoid of rotation with the characteristic small 
semiaxis ^ am — 03/7 and large semiaxis I 3> a™ 
(/ ^ cb/ In 7 am for the tightly-bound states) . Accu- 
rate fitting formulae for energies and other characteris- 
tics of a hydrogen atom in magnetic fields are given in 
[186]. 

The properties of molecules and even the very exis- 
tence of some of their types in strong magnetic fields are 
poorly known, although they have been discussed for al- 
most 40 years. Those diatomic molecules are fairly well 
studied whose axis coincides with the direction of the 
magnetic field. For obvious reasons, the Il2molecule has 
been thoroughly investigated. The approximate formu- 
las for its binding energy at 7 > lO'^ increasing at the 
same rate cx (In 7 )^ as the binding energy of H atoms are 
presented in II6OII . Interestingly, however, numerical cal- 
culations in |187l | , show that this molecule is unstable in a 
moderate magnetic field (in the range 0.18 < 7 < 12.3). 

Al so, t he H2"'"ion has been studied fairly well (see, 
e.g., |188l l): HeH++, H3++, and other exotic one-electron 
molecular ions becoming stable in strong magnetic fields 
were also considered [189j] . 

A strong magnetic field can stabilize polymer molec- 
ular chains aligned along magnetic field lines. These 
chains can then attract one another via dipole-dipole 
interactions and make up a condensed medium. Such 
a po ssibility was first conjectured by Ruderman in 1971 
|185{ . Investigations in the 1980s - 2000s showed that in 
the fields B ~ 10^^-10^'^ G these chains are formed not 
of any chemical elements, but only of the atoms from 
H to C, and undergo polymerization into a condensed 
phase either in a superstrong field or at a relatively low 
temperature, with the sublimation energy of such con- 
dens ate b eing much smaller than predicted by Ruderman 
(see |190j and the references therein). 

The overwhelming majority of researchers of atoms 
and molecules in strong magnetic fields have considered 
them to be at rest. Moreover, in studies of electron 
shells, the atomic nuclei were almost universally assumed 
to be infinitely massive (fixed in space). Such an ap- 
proximation is a gross simplification for magnetic atmo- 
spheres. Astrophysical simulations must take the finite 
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temperature and therefore thermal motion into account. 
Atomic motion across magnetic field lines breaks the ax- 
ial symmetry of a quantum mechanical system. At 7 3> 1 
, specific effects associated with the collective motion of 
a system of charged particles become significant. Specif- 
ically, the decentered states can be populated in which 
electrons are mainly located in a "magnetic well" far 
from the Coulomb center. Th ese e xotic states were pre- 
dicted for hydrogen atoms in [191j |. In the same paper, 
and later in ^192i | , the first studies of the energy spectrum 
of these states were performed. 

Even at low temperatures, when the thermal motion 
of atoms can be neglected in the first approximation, 
the finite mass of the atomic nucleus should be taken 
into consideration, even in a sufficiently strong field; the 
nucleus undergoes oscillations in a magnetic field due to 
Landau quantization even if the generalized momentum 
[l83] describing the motion of the center of masses across 
the field is zero. Different quantum numbers of an atom 
correspond to different vibrational energies that are mul- 
tiples of the cyclotron energy of the atomic nucleus. In 
a superstrong field, this energy becomes comparable to 
the electron shell energies and cannot therefore be dis- 
regarded. 

A comprehensive calculation of hydrogen atom energy 
spectra taking account of motion effects acr oss the strong 
magnetic field was carried out in [193l Il94| , and the cal- 
culation of the probability of different types of radia- 
tive transitions and absorption coefficients in neutron 
star atmospheres in a series of studies was reported in 
[95j . Based on these data, a model of the hydrogen at- 
mosphere of a neutron star with a strong magnetic field 
[162| was elaborated. The database for ast rophy sical cal- 
culations was created using this model in |163| . 

The quantum me chan i cal e ffects of He-|- ion motion 
were considered in [195l . Il96l |. This case is essentially 
different from that of a neutral atom in that the value s 
of the ion generalized momentum are quantized [183j . 
For many-electron atoms, molecules, and ions, the ef- 
fects of motion across the magnetic field remain unex- 
plored. The perturbation the ory applic able to the case 
of small generalized momenta [197l Il98l | may prove suffi- 
cient to si mulate relatively cold atmospheres of neutron 
stars [T6l |. 

6.4 Electron heat and charge transport 
coefficients 

The magnetic field affects the kinetic properties of the 
plasma in a variety of ways (see, e.g., [199], for a review). 
Any magnetic field makes the transfer of charged parti- 
cles (in our case, electrons) anisotropic. It hampers their 
motion and thereby heat and charge transfer by electrons 
in the direction perpendicular to the field, thus generat- 
ing Hall currents. These effects are essential when the 
cyclotron frequency uJc is much higher than the effective 
collision frequency, while the latter remains unaltered in 
a nonquantizing field. 

A quantizing magnetic field exerts a more pronounced 
infiuence on the transfer process. In a weakly quantizing 
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Figure 6: Longitudinal (|]) and transverse (_L) thermal 
conductivities in the iron outer env elope of a neutron 
star at T = 10^ K and B = lO" G [20ll. Solid curves 
are calculations according to |20Clj . dashed lines are the 
classical model, the dotted line is the results neglecting 
thermal averaging. Electrons are degenerate to the right 
of the vertical dashed-dotted line, whose position corre- 
sponds to the equality Tp — T. 



field (in the presence of degeneracy) , kinetic coefficients 
oscillate with variations of matter density about the val- 
ues they would have in the absence of quantization. In a 
strongly quantizing magnetic field, the values of the ki- 
netic coefficients are substantially different from classical 
ones. 



The effects of quantizing magnetic fields on electron 
transfer in plasma have been studied by d iffere nt au- 
thors for many decades (see the references in [199l |). The 
formulas for electron kinetic coefficients of a completely 
ionized plasma convenient to use in astrophysics at ar- 
bitrary p, B, and T were derived in 200]. They were 
used to calculate thermal evolution of neutron stars as 
described in the next section. 

Figure [6] exemplifies the dependences of heat conduc- 
tivity coefficients along and across a magnetic field at the 
plasma characteristics inherent in the blanketing shell of 
a neutron star with field B — 10^'* G; the accurately com- 
puted characteristics (solid curves) are compared with 
the simplified models used in astrophysics previously. 
The dashed lines represent models disregarding Landau 
quantization under the assumption of a strong electron 
degeneracy (to the right of the vertical dashed-dotted 
line) or nondegeneracy (to the left of the vertical line). 
The dotted line is the assumption under which thermal 
scattering of electron energies near the Fermi energy is 
neglected (see ,201.] for the details and references). 
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7 Cooling and thermal radiation 
7.1 Cooling stages 

About 20 s after its birth, a neutron star becomes trans- 
parent to neutrino emission [4^ . carrying away the en- 
ergy to outer space and cooHng the star. Soon after that, 
the temperature distribution in the stellar core charac- 
terized by high heat conductivity reaches equilibrium, 
preserved thereafter throughout the star lifetime (prob- 
ably except for short periods after catastrophic phase 
transitions in the core postulated by certain hypothetical 
models). In line with GR, the equilibrium temperature 
increases toward the center of the star in proportion to 
e~*, where $ is the metric function defined by the star 
hydrostatic model and related to the time-component of 
the metric tensor .Qoq = e ^^, which decreases from the 
surface to the center |202| . 

The stellar crust is for some time hotter than the 
core. The cooling wave reaches the surface within 10- 
100 years; thereafter, the star cools down in the qua- 
sistationary regime where the temperature distribution 
in the heat-insulating blanket at each time instant un- 
ambiguously depends on the core temperature. We note 
that all currently observed neutron stars are at least sev- 
eral centuries old. This means that they are in the state 
of quasistationary cooling in the absence of fast energy 
release in the envelopes. The quasistationarity may be 
disturbed by the explosive thermonuclear burning of ac- 
creted matter [53| or the liberation of energy in the crust 
during starquakes [SS'-'SG^. 

Cooling in the quasistationary regime goes through 
the following stages [SSj . 

1. The neutrino cooling stage lasts ~ 10^ years. Dur- 
ing this time, the core cools largely via neutrino 
emission in various physical reactions (30| , the main 
ones being direct (if present) and modified Urea pro- 
cesses (depending on the particles involved), as well 
as neutrino bremsstrahlung radiation. 

2. The photon cooling stage is the final one. It begins 
at the stellar age t > 10^ years when the lowered 
core temperature makes neutrino emission (strongly 
temperature-dependent) weaker than in cooling via 
heat transfer through the envelope and conversion 
into surface electromagnetic radiation. 

The cooling curve of a neutron star depends on its 
mass M; the model of superdense matter in the core 
determining the equation of state (hence, the radius R) 
and composition of the core (hence, the intensity of neu- 
trino emission at a given mass); and the envelope proper- 
ties: (a) thermal conductivity determining at a given 
core temperature, (b) neutrino luminosity in the crust, 
(c) sources of heating and their intensity. Characteristic 
thermal conductivity and neutrino luminosity of the en- 
velopes at each time instant t (i.e., at the model-specific 
temperature T distribution in the envelopes), in turn, 
depend on the stellar mass M, the radius R, and the 
the magnetic field (both magnetic induction B and the 
configuration of magnetic force lines may be essential). 



Comparing the observed Lj and t for neutron stars 
with the cooling curves allows estimating AI and R 
and placing bounds on the theoretical models of super- 
dense matter. This method for parameter evaluation is 
largely applicable to isolated neutron stars. By contrast, 
most neutron stars in binary systems have an additional 
source of energy (accretion) and an additional source of 
X-ray radiation (accretion disk), often much more pow- 
erful than Lj. 

7.2 Thermal structure 

The complete set of equations describing the mechanical 
and thermal structure and evolution of a spherically sym- 
metric star at hydrostatic equilibrium in the framework 
of GR was obtained by Thorne j202| . These equations 
are easy to transform to the form holding for the stellar 
envelope with radial heat transfer, a smooth temperature 
distribution over the surface, and a force-free magnetic 
field. Under the assumption that heat transfer and neu- 
trino emission are quasistationary, these equations re- 
duce to a system of ordinary differential equations for the 
metric function $, local density of the radial heat flow 
Fr, temperature T and gravitational mass m comprised 
inside a sphere of radius r, as functions of pressure P 
(e.g., ;204]). GR correction factors in this system depend 
on the mass fraction (M — m) /Moutside the equipoten- 
tial surface being considered and on the P/ (pc^) ratio en- 
tering the TOV equation. At the outer-inner crust inter- 
face, {M-m)/M - 10-5 and P/{pc^) ^ 10"^, therefore, 
the GR correction factors in the outer shells are almost 
constant. Bearing this constancy in mind and disregard- 
ing the geometric thickness of the heat-insulating layer 
compared with _R, in the absence of heat sources and 
sinks in the blanketing envelope of a neutron star, we 
obtain that the radial heat flow Fr is constant and equal 
to ctsbT's'j where ass is the Stefan-Boltzmann constant. 
Here and hereafter, we discriminate between the local 
surface temperature and integral effective tempera- 
ture Tcff, because Tg may vary over the surface. Under 
the above conditions, calculating the ther mal s tructure 
amounts to solving a simplified equation |205| . which 
can be written in the same form as the nonrelativistic 
equation ndT/dP — Fr/{pg). Such an approximation is 
used in the majority of neutron star cooling research [88| . 
However, since magnetars have stronger magnetic fields 
and surface luminosities than ordinary neutron stars and 
have, in addition, internal sources of energy, one has to 
solve for them the complete set of equations, taking neu- 
trino emission rate per unit volume Qi, and heat sources 
Qh into account, instead of the simplified heat transfer 
equation .100 . 204] . 

The effective radial thermal conductivity at a local 
surface area in a magnetic field is k = K|| cos^ 9b + 
K± sin^ 9b, where 9b is the angle between magnetic force 
lines and the normal to the surface, and K|| and kj_ are 
components of the heat conductivity tensor responsible 
for the transfer along and across the force lines. In the 
heat-insulating envelope of a neutron star each of the 
components /t|| and k± contains radiative Kr and elec- 
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tron Ke constituents. Photon heat conduction prevails 
{Kr > Ke) in the outcrmost (typically nondegenerate) 
layers, while electron heat conduction plays the main 
role in deeper, moderately or strongly degenerate layers. 
The total thermal flux along a given radius r (the local 
luminosity related to thermal but not neutrino losses) 
is defined by the flux density integral over the sphere 
of this radius, Lr = J sinOdO dip r^Fr{0,ip), where {0,(p) 
are respectively the polar and azimuthal angles. 

The use of equations holding for a spherically symmet- 
ric body at each point of the surface assumes that the 
mean radial temperature gradient is much greater than 
the lateral one. The estimates made in [204| indicate that 
this condition is fulfilled with a good accuracy at the 
largest part of the star surface, and corrections for de- 
viations from a one-dimensional approximation make a 
negligibly small contribution to the total luminosity; this 
allows disregarding them in the first approximation. 

In the quasistationary regime, the temperature of a 
neutron star increases monotonically from the external 
layers of the atmosphere to the interior of the envelope 
until it reaches equilibrium (usually in the outer crust). 
However, magnetars must have sources of heating in the 
envelopes capable of maintaining their high luminosity; 
for this reason, temperat ures profiles in magnetar en- 
velopes are nonmonotonic [lOOj . 

7.3 Cooling curves 

The nonstationary proble m is described by the same 
thermal balance equations |204| . but the difference Q — 
Qu — Qh is supplemented by the term Ce~* dT/dt, where 
C is the heat capacity per unit volume '203]. Strictly 
speaking, Q should be additionally supplemented by yet 
another term describing the release of latent melting heat 
during the movement of the Coulomb fluid-crystal inter- 
face with a change in temperature. But this term is al- 
ways neglected in the available programs for calculating 
the thermal evolution of neutron stars. Following the 
classical work [2051] , the nonstationary problem is solved 
in the interior of a neutron star where density surpasses 
a certain threshold value pb, while for external envelopes 
at p < Pb, whose relaxation time is short relative to char- 
acteristic times of thermal evolution, a stationa ry sy stem 
of equations is solved. Traditionally following [205| . one 
chooses Pb — 10^° g cm~^, but sometimes other pb values 
prove more suitable, dep ending on the concrete problem 
of interest [100, 204|2Ql[2^. The relation between the 
heat flow across the boundary pb and the temperature 
Tb at this boundary obtained by solving the stationary 
problem for the envelopes serves as a boundary condition 
for the nonstationary problem in the internal region. It 
primarily depends on heat conductivity in the sensitivity 
strip on the p — T plane near the turning point,''^ where 
Ke ~ Kr |205l ] . Analytic estim ates for the position of this 
point were obtained in [201| . 

The quantities Tg, Teff, and are defined in the local 
reference frame at the neutron star surface. The "appar- 
ent" quantities in the frame of a distant observer should 
be corrected for the redshift (see Sect. 12. ip . The solution 
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Figure 7: Cooling curves of neutron stars compared with 
some obs ervat ion-based estimates of their temperatures 
and ages [208j . The cooling curves for different models 
of chemical composi tion of the heat-insulating envelope 
(in accordance with '2061]) correspond to different accu- 
mulated masses AM: of light elements: solid and dotted 
curves correspond to an iron envelope and an envelope of 
lighter nuclear composition, dashed curves correspond to 
a partly substituted envelope. The upper three curves 
correspond to a star with the mass M = 1.3 Mq, un- 
dergoing standard cooling by Murca processes and the 
lower three, to a star with the mass M —1.5 AIq, under- 
going enhanced cooling by direct Urea processes. Dots 
with error bars correspond to the published estimates of 
ages t and effective surface temperatures T^ of neutron 
stars; arrows pointing down indicate the upper bounds 
onT,^. 



of th e cooling problem is described at greater depth in 
|203j (see also reviews [H, l208j | and references therein) . 

The envelopes of a neutron star at birth consist of iron- 
group elements, which explains why calculations of cool- 
ing were for a long time made for iron-rich shells alone. 
However, the envelopes of a star that has passed through 
an accretion stage may consist of lighter elements. Ac- 
creted envelopes have a higher electron conductivity 
than iron-rich ones because more weakly charged ions 
less effectively scatter electrons. In other words, ac- 
cretion makes the envelopes more "transparent" to the 
passing heat 209]. The core temperature at the neu- 
trino cooling stage is regulated by neutrino emission and 
is practically independent of the properties of the en- 
velopes; therefore, their transparency makes the star 
brighter due to an enhanced T^g. At the later pho- 
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Figure 8: The same as in Fig. [71 but for the iron envelope 
and different magnetic fields: solid and dotted curves 
correspond to B = and 10^^ G respectively, while the 
dashed curves correspond to an intermediate case. 



ton stage, the transparent envelopes more readily trans- 
mit heat and the star fades away faster. These effects 
are especially demonstrative when comparing the solid, 
dashed, and dotted curvesQ in Fig. [71 the difference be- 
tween which is due to the different mass of accreted mat- 
ter AM. 

Similarly, cooling depends on a superstrong magnetic 
field. In a strong magnetic field where the electron cy- 
clotron frequency w exceeds the typical frequency of 
their collisions with plasma ions, heat transfer across 
magnetic field lines is hampered; therefore, those regions 
in which the lines are directed toward the surface be- 
come cooler. Oscillations of the heat conductivity co- 
efficients (see Fig. [B]), caused by Landau quantization 
facilitate heat transfer along magnetic force lines on an 
average, making the regions near the magnetic poles hot- 
ter. Taken together, enhanced luminosity near the poles 
and its reduction at the equator make integral luminos- 
ity of a star in a moderate magnetic field B < 10^^ G 
virtually the same as in the absence of a magnetic field. 
Although the temperature distribution over the surface 
depends on the field strength and configuration, the in- 
tegral luminosity is virtually unrelated to B for a moder- 
ate dipole field |206] as well as for a moderate small-scale 
magnetic field [2 111] . However, in the superstrong field of 



•^The cooling curves pre sente d in Figs. [7] and [8] were calculated 
by D G Yakovlev for Ref. I210II using a relatively soft equation of 
state , fo r which the direct Urea processes open at M > 1.462 Mq 
(see m). 



magnetars, B > 10^'* G, the enhancement of the trans- 
parency near the poles is more significant, which causes 
the effect of the integral transparency enhancement of 
the envelopes, analogous to the effect of accretion, as 
shown in Fig. [H 

7.4 Effective temperatures 

The upper and lower groups of three curves in each of 
Figs. [7] and [8] correspond to standard and enhanced cool- 
ing. The latter occurs when direct Urea processes oper- 
ate at the neutrino cooling stage in a star of a sufficiently 
large mass. The dots with error bars indicate estimated 
ages t and effective surface temperatures Tjg obt ained 
from observational data summarized in Ref. j208l ]. We 
see that under favorable conditions, cooling curves give 
an idea of the stellar mass and properties of envelopes: 
the coldest stars of a given age appear to undergo en- 
hanced cooling and are therefore massive, whereas the 
hottest ones have accreted envelopes. 

Enhanced cooling may be a consequence not only of 
direct Urea processes in npe/i-matter but also of analo- 
gous hyperon and quark Urea process es in exotic models 
of the inner core of a neutron star [125| . The rate of 
direct Urea processes is limited by a gap in the energy 
spectrum of superfluid nucleons [20 ] ; therefore, nucleon 
superfiuidity smooths the dependence of cooling curves 
on stellar mass, making the "weighing" of the st ar by 
measuring its effective temperature more feasible [208j . 
Moreover, nucleon superfluidity, as well as hyperon and 
quark (color) superfiuidity in exotic models, decreases 
the heat capacity of the stellar core [2^, which also af- 
fects cooling 208]. 

Thus, comparison of measured ages and temperatures 
with cooling curves allows one, in principle, to deter- 
mine stellar mass, or to conjecture the composition of 
the core and heat-insulating envelopes when the mass 
is known from independent estimates. But the results 
thus obtained should be interpreted with caution for the 
following reasons. Effective temperatures are typi- 
cally measured by varying the parameters of a theoretical 
model used to calculate the emission spectrum. The pa- 
rameters are chosen so as to most accurately reproduce 
the observed spectrum, but the final result strongly de- 
pends on the choice of the model. Certain estimates 
Tjjg, presented in Figs. [71 and El were obtained using 
models of nonmagnetic and magnetic hydrogen atmo- 
spheres, while others simply assume that radiation is 
described by the Planck spectrum. For example, the 
fit to the three-component model spectrum was used for 
the thermal component of the PSR BI055-52 spectrum: 
the power-law component was added to two black-body 
components, whose cooling is believ ed to be responsible 
for heat radiation from the surface [212l ]. The result of 
this fitting is marked in the figures by the numeral 1, 
while the one marked by 2 was obtained for the radio- 
quiet neutron star RX J1856. 4-375 4 bas ed on a physical 
model of t he m agnetic atmosphere |213| (we will discuss 
it in Sect. [7T5] in more detail). For comparison, the au- 
thors of [213»] also fitted the observed X-ray spectrum 



PHYSICS OF NEUTRON STARS 



23 



of RX J1856.4-3754 with the Planck spectrum. When 
plotted in Figs. [7] and [3 the result of such a simplis- 
tic fitting would almost coincide with the point 1. This 
means that the systematic error in the position of point 
1 resulting from the absence of a model for the forma- 
tion of the PSR B1055-52 spectrum may be of the same 
order as the distance between points 1 and 2, that is, 
significantly greater than the statistical fitting error. 

7.5 Masses and radii 

Analysis of the thermal radiation spectrum of a neutron 
star provides information not only on its effective tem- 
perature but also on the radius and mass. Let us first 
consider the Planck spectrum disregarding interstellar 
absorption and possible non uniformity of the tempera- 
ture distribution over the surface. The position of the 
spectral maximum gives T^, and the total (bolometric) 
flux of incoming radiation is found from its mea- 
sured intensity. If a star is at a distance D, its apparent 
photon luminosity is obtained as = AnD^Fy^^. On 
the other hand, according tot he Stefan-Boltzmann law 
= Attosb Rto C^^)"^ ^ which allows one to estimate 

In reality, the comparison of theoretical and measured 
spectra actually involves more unknown variables. The 
spectrum is distorted by absorption in the interstellar 
gas; therefore, spectral analysis can be used to deter- 
mine average gas concentration over the line of sight. 
If the distance D is unknown, it can be estimated un- 
der the assumption of the typical concentration of the 
interstellar gas in a given galactic region, using D as 
a fitting parameter. The temperature distribution over 
the stellar surface may be nonuniform. For example, 
the thermal spectrum in the presence of hot polar caps 
consists of two blackbody components, each having its 
own values of T^^ and Roo- Finally, because the star 
is not a perfect black body, the real radiation spectrum 
differs from the Planck spectrum. Spectrum simulation 
is a difficult task, involving solution of the equations of 
hydrostatic e quilib rium, energy balance, and radiation 
transmission |214| . Coefficients of these equations de- 
pend on the chemical composition of the atmosphere, 
effective temperature, acceleration of gravity, and mag- 
netic field. Different assumptions of the chemical com- 
position, Tcff, Zg, and B, values lead to different model 
spectra; their comparison with the observed spectrum 
permits obtaining acceptable parameter values. Know- 
ing the shape of the spectrum, allows calculating the 
fitting coefficient by the Stefan-^Boltzmann formula and 
finding Roc from i^j^j. The simultaneous finding of Zg 
and Roo = ^^(1 + Zg) allows one to calculate the mass Af 
on the base of Eqs. ([T]) and ([3]). 

Let us consider some of the problems arising from 
the estimation of parameters of neutron stars from their 
observed thermal spectra using the "Walter star" RX 
J1856. 4-3754 as an example. It is a nearby radio-quiet 
neutron star, discovered in 1996 as a soft X-ray s ourc e 
[OOj and identified a year later in the optical range [215{ . 
The first parallax measurement by the "planetary cam- 



era" (PC) on board the H ubble space observatory gave 
the value Z? w 60 pc [216j , which corresponded to a very 
small radius R. A later improvement of the distance to- 
gether with a preliminary spectral analysis gave a some- 
what larger R which might correspond to a quark star 
|217||. A subsequent reanalysis of the data gave D « 120 
pc '218] , and independent treatment of the same data by 
different authors has led to D w 140 pc ^219]. Finally, 
the measurement of parallax by the high resolution cam- 
era (HRC) of the same observatory yielded D « 160 pc. 
At the same time, it turned out that the spectrum of the 
Walter star is not described by the blackbody model: the 
fit of its X-ray region to the Planck spectrum predicts 
a much weaker luminosity in the optical range than the 
observed one. Attempts to describe the observed spec- 
trum in terms of the models of atmospheres of different 
chemical compositions without a magnetic field a nd of 
the two-component Planck model are reported in [22 Ij 
and [222], respectively. It turned out that the hydrogen 
atmosphere model reproducing the X-ray spectral region 
predicts too high a luminosity in the optical range and 
models of atmospheres of a different chemical compo- 
sition predict absorption lines unseen in observational 
data. Fitting to the two-component model for the soft 
component of the composite spectrum leads to a bound 
on the radius Roo > 17 km x(D/120 pc), which is diffi- 
cult to relate to theoretical calculations of neutron star 
radii. 

A simulation of the neutron star spectrum based on 
the solution of a system of equations for radiation trans- 
fer in a partly ionized hydrogen atmosphere of finite 
thickness above the condense d sur face in a strong mag- 
netic field was proposed in 
the atmosphere model from 



215 1 . The authors used 
162| . based on the equa- 



tion of state of the hydrogen plasma in a strong mag- 
netic field and absorption/scattering coefficients in such 

a plasma presented in [95] . At S (3-4) x 10^^ G, 
TOO _ /-/I Q/i _i_ n n^A V 1 no Ik' r, —no 



= (4.34 ± 0.03) X W K, Zg = 0.25 ± 0.05, and 

-Roo — 17.2^01 duo km, they managed for the first time 
to reproduce the measured spectrum of RX J1856. 4-3754 
in the X-ray to optical frequency range within the mea- 
surement errors of the best space and terrestrial obser- 
vatories. Here, the errors are given at the significance 
levels Icr, and duo ^ -D/(140 pc). Taking relations 
11]) -([SI) into account, one founds from these estimates, 
that for this neutron star R — 13.8^0 e rfi4o km and 



M — 1.68_o i5 di4o -^0- Forgetting for a moment the 
multiplier c?i4o, one might conclude that the 68% confi- 
dence area lies above all theoretical dependences R{M), 
shown in Fig. H] It could mean that either the measure- 
ments or the theoretical model are not accurate (if not 
to consider a possibility of a superstiff equation of s tate, 
not shown in Fig. El). The estimate D sa 160 pc [22oj 
shifts the values of R and M still farther from the the- 
oretical dependences R{M). Moreover, such a massive 
star should have undergone enhanced cooling, which is 
not observed in Figs. [7] and [HI However, confirmati on of 
the estimate D « 120 pc [218] in a recent paper [223j 
eliminates these contradictions. 
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Thus, the uncertainty of distance measurements 
proves more important than the inaccuracy of spectrum 
fitting. An even greater uncertainty is associated with 
the choice of a theoretical model. For example, fitting 
the X-ray region of the RX J1856. 4-375 4 spe ctrum with 
the black-body spectrum, presented in |213l l for means 
of comparison, gives Roo ~ 5 duo km. 

Similar problems are encountered in the analysis of all 
known thermal spectra of isolated neutron stars. They 
are not infrequently supplemented by uncertainties of 
spectrum division into thermal and nonthermal compo- 
nents (see 89] and the references therein). 



8 Conclusions 



Neutron stars are miraculous objects in which Nature 
assembled its puzzles, whose solution is sought by seem- 
ingly unrelated branches of science, such as the physics 
of outer space and the micro-world, giant gravitating 
masses, and particle particles. This makes neutron stars 
unique cosmic laboratories for the verification of basic 
physical concepts. In the last 50 years, both theoretical 
and observational studies of neutron stars have been de- 
veloping at a progressively faster pace following advances 
in nuclear and elementary physics on the one hand, and 
astronomy and experimental physics on the other hand. 

The present review outlines some aspects of neutron 
star physics, describes methods for measuring their tem- 
perature, masses, and radii, and illustrates the relation 
between the theoretical interpretation of these data and 
the solution of fundamental physical problems. 

The work was supported by the Russian Foundation 
for Basic Research (grant 08-02-00837) and the State 
Program for Support of Leading Scientific Schools of 
the Russian Federation (grant NSh-3769.2010.2). 



Note added in proof 

When this paper was being prepared for publication, P 
Demorest et al. reported the record-breaking mass M = 
1.97 ± 0.04, Mq of the neutron star in the PSR J1614- 
2230 binary system [224]. This estimate was obtained by 
measuring the Shapiro delay parameters (see Section 2.1 
of the present review). Plotting it in our Fig. 4 results in 
appearance of the corresponding vertical strip slightly to 
the left of point 2. In conformity with the discussion in 
Section 4.4, it prompts that superdense matter cannot 
be characterized by equations of state softer than Sly. 
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